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[i]  The  LaRC  Lagrangian  Chemical  Transport  Model  (LaRC  LCTM)  is  used  to  simulate 
the  kinematic  and  chemical  evolution  of  an  ensemble  of  trajectories  initialized  from 
Halogen  Occultation  Experiment  (HALOE)  and  Polar  Ozone  and  Aerosol  Measurement 
(POAM)  III  atmospheric  soundings  over  the  SAGE  III — Ozone  Loss  and  Validation 
Experiment  (SOLVE)  campaign  period.  Initial  mixing  ratios  of  species  which  are  not 
measured  by  HALOE  or  POAM  III  are  specified  using  sunrise  and  sunset  constituent  CH4 
and  constituent  PV  regressions  obtained  from  the  LaRC  IMPACT  model,  a  global  three 
dimensional  general  circulation  and  photochemical  model.  Ensemble  averaging  of  the 
trajectory  chemical  characteristics  provides  a  vortex-average  perspective  of  the 
photochemical  state  of  the  Arctic  vortex.  The  vortex-averaged  evolution  of  ozone,  chlorine, 
nitrogen  species,  and  ozone  photochemical  loss  rates  is  presented.  Enhanced  chlorine 
catalyzed  ozone  loss  begins  in  mid-January  above  500  K,  and  the  altitude  of  the  peak  loss 
gradually  descends  during  the  rest  of  the  simulation.  Peak  vortex  averaged  loss  rates  of  over 
60  ppbv/day  occur  in  early  March  at  450  K.  Vortex  averaged  loss  rates  decline  after  mid- 
March.  The  accumulated  photochemical  ozone  loss  during  the  period  from  1  December 
1999  to  30  March  2000  peaks  at  450  K  with  net  losses  of  near  2.2  ppmv.  The  predicted 
distributions  of  CH4,  03,  denitrification,  and  chlorine  activation  are  compared  to  the 
distributions  obtained  from  in  situ  measurements  to  evaluate  the  accuracy  of  the 
simulations.  The  comparisons  show  best  agreement  when  diffusive  tendencies  are  included 
in  the  model  calculations,  highlighting  the  importance  of  this  process  in  the  Arctic  vortex. 

Sensitivity  tests  examining  the  large-scale  influence  of  orographically  generated  gravity- 
wave  temperature  anomalies  are  also  presented.  Results  from  this  sensitivity  study  show 
that  mountain-wave  temperature  perturbations  contribute  an  additional  2-8%  03  loss 
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1.  Introduction 

[2]  Arctic  springtime  ozone  loss  has  been  observed  for  a 
number  of  years  [e.g.,  McKenna  et  al ,  1990;  Salawitch  et 
al ,  1990;  Hoffman  and  Deshler ,  1991;  Proffitt  et  al ,  1993] 
although  the  total  reduction  in  column  ozone  amounts  has 
typically  been  much  less  than  for  the  Antarctic  region.  The 
absence  of  extensive  Northern  Hemisphere  ozone  depletion 
is  primarily  due  to  the  more  transient  nature  of  the  Arctic 
vortex,  which  tends  to  experience  some  degree  of  disruption 
over  the  course  of  the  winter,  leading  to  warmer  temper¬ 
atures,  and  less  persistent  polar  stratospheric  cloud  (PSC) 
decks.  The  lack  of  persistent  PSCs  means  that  the  denitri- 
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fication  necessary  for  maintaining  high  levels  of  chlorine 
activation  required  for  significant  spring  time  chlorine- 
catalyzed  ozone  destruction  is  often  not  realized  in  the 
Arctic  l Solomon ,  1999].  However,  long-term  observations 
of  total  ozone  by  the  Total  Ozone  Mapping  Spectrometer 
(TOMS)  indicate  significant  declines  in  total  ozone  over  the 
Arctic  during  the  last  decade  [Newman  et  al.,  1997]. 
Column  ozone  losses  in  the  lower  stratosphere  of  as  much 
as  50%  of  the  column  above  the  150  hPa  level  have  been 
calculated  by  Muller  et  al  [1996]  using  coincident  measure¬ 
ments  of  03  and  CH4  by  HALOE.  The  potential  severity  of 
Arctic  ozone  loss  was  made  quite  evident  during  1997, 
when  an  anomalously  persistent  vortex  resulted  in  record 
low  ozone.  In  contrast,  the  1998  boreal  winter  showed 
warmer  vortex  temperatures  and  more  normal  springtime 
column  ozone  amounts  [Schultz  et  al .,  2001]. 

[3]  The  wide  variability  in  Northern  Hemisphere  spring¬ 
time  ozone  amounts  is  due  to  a  strong  coupling  between 
chlorine  loading  and  planetary-wave  dynamics.  Interannual 
variability  in  the  Northern  Hemisphere  circulation  leads  to 
significant  long-term  variability  in  Arctic  temperatures, 
ozone  transport  and  photochemistry.  Episodic  planetary- 
wave  disturbances  lead  to  significant  intraseasonal  ozone 
transport  and  photochemical  variability  within  the  Arctic. 
The  transient  nature  of  the  Arctic  vortex  also  means  that 
variations  in  vertical  transport  and  diffusive  mixing  pro¬ 
cesses  tend  to  play  a  more  important  role  in  determining  the 
distribution  of  column  ozone  in  the  Arctic  than  in  the 
Antarctic,  where  the  vortex  is  more  isolated  and  stable. 
This  transience,  coupled  with  wide  variations  in  the  extent 
of  denitrification  and  chlorine  activation,  tends  to  compli¬ 
cate  the  observationally  based  determination  of  the  chemical 
loss  within  the  Arctic  vortex. 

[4]  The  1999/2000  winter  presents  a  unique  opportunity 
to  test  our  understanding  of  Arctic  ozone  loss  because  the 
extent  and  duration  of  the  Arctic  cold  pool  was  one  of  the 
largest,  and  longest  on  record  [Manney  and  Sabutis , 
2000]  therefore  significant  denitrification  and  chlorine 
activation  should  have  occurred.  Furthermore,  the  large 
amount  of  in  situ,  remote,  and  balloon  data  available  from 
the  SAGE  III  Ozone  Loss  and  Validation  Experiment 
(SOLVE)  Campaign,  which  was  conducted  during 
November- March  1999/2000,  provides  an  unprecedented 
opportunity  to  test  our  theoretical  understanding  of  Arctic 
ozone  loss  processes. 

[5]  The  LaRC  Lagrangian  Chemical  Transport  Model 
(LaRC  LCTM)  [Pierce  et  al ,  1996,  1999]  was  used  during 
the  1999/2000  SOLVE  campaign  to  provide  real-time  fore¬ 
casts  of  the  hemispheric  distribution  of  catalytic  ozone  loss 
processes,  transport,  and  mixing  in  the  lower  and  middle 
stratosphere.  The  model  uses  daily  satellite  observations  of 
ozone  and  other  key  species  from  the  Halogen  Occultation 
Experiment  (HALOE)  [Russell  et  al ,  1993]  and  the  Polar 
Ozone  and  Aerosol  Measurement  (POAM  III)  [Lucke  et  al , 
1999]  to  initialize  an  ensemble  of  trajectories.  The  LaRC 
LCTM  then  predicts  the  photochemical,  microphysical,  and 
diffusive  tendencies  along  these  trajectories.  The  number  of 
trajectories  increases  as  the  integration  proceeds,  resulting 
in  an  ensemble  Lagrangian  prediction  of  the  chemical  and 
microphysical  state  of  the  Arctic  stratosphere.  The  SOLVE 
LaRC  LCTM  predictions  provide  a  large-scale  perspective 
of  ozone  loss  and  transport  processes  to  aid  in  the  inter¬ 


pretation  of  in  situ  and  remote  data  collected  during  the 
SOLVE  field  campaign.  Satellite  solar  occultation  measure¬ 
ments,  in  conjunction  with  the  LaRC  LCTM,  significantly 
enhance  the  vertical  and  longitudinal  extent  of  the  observa¬ 
tional  and  theoretical  database  for  SOLVE  and  allow  us  to 
examine  the  relative  contributions  of  advective,  diffusive, 
and  photochemical  processes  in  determining  the  observed 
evolution  of  ozone  within  the  Arctic  vortex. 

[6]  The  format  of  this  paper  is  as  follows.  We  begin  with 
a  discussion  of  recent  modifications  to  the  model  which 
were  needed  to  properly  simulate  the  Arctic  stratosphere 
and  a  description  of  the  SOLVE  simulation.  Then  we 
present  a  discussion  of  the  vortex-averaged  evolution  of 
the  simulation  over  the  course  of  the  1999/2000  winter. 
Next  we  present  comparisons  of  the  model  results  with  in 
situ  observations  from  instruments  onboard  the  ER2,  exam¬ 
ine  two  specific  case  studies  in  detail,  and  summarize  the 
integrated  budgets  of  03  and  NOy  within  the  Arctic  vortex. 
Finally,  we  present  a  sensitivity  study  that  explores  the 
impact  of  orographically  generated  gravity-wave  temper¬ 
ature  perturbations  on  the  vortex  03  and  NOy  budgets. 

2.  The  LaRC  Lagrangian  Chemical  Transport 
Model  (LaRC  LCTM) 

[7]  The  LaRC  LCTM  is  formulated  in  isentropic  coor¬ 
dinates  and  uses  once  daily  12Z  winds  and  temperatures 
obtained  from  the  United  Kingdom  Meteorological  Office 
(UKMO)  Upper  Atmospheric  Research  Satellite  (UARS) 
assimilation  [Swinbank  and  O’Neill ,  1994].  The  horizontal 
resolution  of  the  UKMO  data  is  3.75°  in  longitude  by  2.5° 
in  latitude  and  has  18  pressure  levels  from  the  surface  to 
0.4hPa.  The  vertical  resolution  is  5  pressure  levels  per 
decade  that  are  uniform  in  log  pressure  (approximately 
3km  resolution).  Assimilated  winds  and  temperatures  are 
interpolated  to  34  isentropic  surfaces  extending  from  240K 
to  400K  by  10K,  400K  to  550K  by  25K,  600K  to  1000K  by 
100K,  and  1000K  to  2000K  by  200K  increments.  This  gives 
a  nominal  vertical  resolution  for  the  trajectory  calculations 
of  approximately  1.0  km  in  the  lower  stratosphere. 

[8]  The  UKMO  assimilated  temperatures  and  the  MID¬ 
RAD  radiative  code  [Shine,  1987]  are  used  to  compute 
radiative  heating  rates  (RtotaiX  which  are  used  to  determine 
the  vertical  displacement  of  the  trajectories.  Mass  conser¬ 
vation  is  assured  by  adjusting  the  short-wave  component  of 
the  radiative  heating  (RShortwave)  such  that  the  globally 
averaged  radiative  heating  rate  on  a  pressure  surface  is 
zero.  The  adjustment  (1.  —  Rtotai/Rshortwave)  scales  the 
shortwave  component  of  the  heating  and  leaves  the  long¬ 
wave  component  of  the  heating  unmodified.  This  adjust¬ 
ment  procedure  was  chosen  since  it  tends  to  maintain  the 
strong  radiative  cooling  found  in  the  Arctic  vortex  during 
polar  night  much  better  than  a  global  adjustment  of  the  total 
heating  rate  and  results  in  better  predictions  of  the  vertical 
profiles  of  long-lived  species  with  strong  vertical  gradients 
(such  as  N20  and  CH4)  in  the  Arctic  vortex. 

[9]  The  isentropic  data  base  for  the  model  also  includes 
daily,  3D  distributions  of  winds,  temperatures,  potential 
vorticity,  Q  [Fairlie  et  al ,  1999],  a  measure  of  the  relative 
strength  of  strain  and  rotation  in  the  flow,  stream  function, 
velocity  potential,  overhead  column  ozone,  and  a  circulation 
flag  which  identifies  the  polar  vortex  (regions  where  Q  is 
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negative  and  vorticity  is  positive  when  integrated  along 
streamlines  defined  by  the  rotational  component  of  the 
flow)  as  well  as  the  location  of  closed  anticyclonic  circu¬ 
lations  (regions  where  Q  is  negative  and  vorticity  is  neg¬ 
ative  when  integrated  along  rotational  streamlines).  The 
overhead  column  ozone  is  used  for  the  model  photolysis 
rates  and  is  determined  from  PV-theta  mapped  ozone  dis¬ 
tributions  that  combine  information  from  HALOE,  SAGE 
II,  and  POAM  III  ozone  profiles  observed  during  SOLVE. 
Separate  PV-03  regressions  are  obtained  by  binning  the 
satellite  measurements  in  PV  bins  on  each  isentropic  level 
for  3  distinct  Q  based  circulation  types  (vortex,  anticyclone, 
or  neither)  for  both  Northern  and  Southern  Hemispheres, 
and  3  different  time  periods  (October-November  1999, 
December- January  2000,  and  February -March  2000). 
Comparison  of  the  PV-theta  mapped  column  03  with 
TOMS  shows  good  agreement  in  the  daily  evolution  of 
the  large-scale  features  although  the  PV-theta  mapped 
column  is  lower  than  TOMS  in  regions  with  particularly 
high  column  ozone.  This  is  likely  due  to  the  fact  that  very 
few  of  the  satellite  measurements  go  below  400K-425K 
and  so  the  ozone  column  in  the  lower  most  stratosphere  is 
missed.  Comparisons  with  other  PV-theta  mapped  03  prod¬ 
ucts  [Randall  et  al .,  2002]  shows  very  good  quantitative 
agreement. 

[10]  The  chemical  module  of  the  LaRC  LCTM  is  from  the 
NASA  Langley  Interactive  Modeling  project  for  Atmos¬ 
pheric  Chemistry  and  Transport  (IMPACT)  model  [Al-Saadi 
et  al ,  2002;  Pierce  et  al ,  2000;  Eckman  et  al ,  1995],  which 
is  a  coupled  3D  chemical-dynamical  model  of  the  tropo¬ 
sphere  and  stratosphere.  The  implementation  of  the  chem¬ 
ical  module  in  the  LaRC  LCTM  is  described  by  Pierce  et  al 
[1999]  and  recent  developments  are  described  here. 

[11]  The  representation  of  stratospheric  heterogeneous 
chemical  processes  and  tropospheric  chemistry  have  been 
updated,  allowing  us  to  further  improve  agreement  of 
modeled  ozone  with  observations  in  the  troposphere  and 
lower  stratosphere.  With  these  changes,  35  chemical  fami¬ 
lies  and  individual  species  are  explicitly  transported.  A  total 
of  53  chemical  species  are  followed  using  41  photolytic 
reactions  and  101  gas-phase  reactions.  Concentrations  of 
species  not  transported  are  obtained  by  partitioning  of  the 
chemical  families  using  standard  photochemical  equilibrium 
techniques.  The  parameterization  of  heterogeneous  process¬ 
ing  has  been  improved  by  considering  a  more  comprehen¬ 
sive  reaction  set  and  by  using  the  analytical  scheme  of 
Carslaw  et  al  [1995]  to  determine  the  sulfate  aerosol 
composition  and  effective  reaction  rates.  Nine  heterogene¬ 
ous  reactions,  involving  N205,  chlorine,  and  bromine  spe¬ 
cies,  are  now  considered  on  both  liquid  aerosol  and  polar 
stratospheric  cloud  surfaces.  Denitrification  is  parameter¬ 
ized  as  sedimentation  of  HNO3  in  PSC  particles  in  a  manner 
similar  to  Chipperfield  [1999].  The  amount  of  condensed- 
phase  HN03  is  determined  through  equilibrium  relation¬ 
ships  with  the  presence  of  condensed-phase  HN03  is  taken 
as  an  indicator  of  type  I  PSCs,  with  an  associated  weak 
removal  time  constant  of  25  days;  the  presence  of  con¬ 
densed-phase  H20  is  taken  as  an  indicator  of  type  II  PSCs, 
with  a  much  faster  removal  time  constant  of  2/3  day.  These 
time  constants  are  derived  assuming  mean  fall  velocities  of 
0.04  km/day  and  1.5  km/day  for  the  nominal  1.0  km 
vertical  resolution  of  the  isentropic  database.  These  fall 


velocities  are  appropriate  for  particles  with  1.0 -2.0  and 
20.0  micron  diameters,  respectively  [Toon  et  al,  1989]. 
When  type  II  PSCs  are  present,  the  local  vapor-phase  H20 
concentration  is  reduced  accordingly,  and  dehydration 
through  sedimentation  is  represented  using  the  2/3  day 
removal  time  constant. 

[12]  Accounting  for  an  exchange  of  condensed-phase 
HN03  and  H20  between  trajectories  through  sedimentation 
and  evaporation  presents  some  unique  challenges  within  the 
Lagrangian  framework  of  the  LaRC  LCTM.  Vertical 
exchange  of  condensed  HN03  and  H20  is  evaluated  by 
determining  the  average  amount  of  condensed  HN03  or 
H20  which  has  sedimented  from  trajectories  above  into 
each  model  gridbox  (defined  by  the  horizontal  and  vertical 
resolution  of  our  isentropic  meteorological  fields)  and  then 
adding  this  amount  to  the  gas-phase  HN03  and  H20 
concentrations  of  all  trajectories  within  that  grid  box.  New 
condensed-phase  HN03  and  H20  concentrations  are  then 
evaluated  for  each  trajectory  within  the  grid  box  using 
equilibrium  relationships,  followed  by  sedimentation  of 
any  new  condensed  HN03  or  H20  into  the  next  vertical 
grid  box.  This  sequential  procedure  accounts  for  sedimen¬ 
tation  through  multiple  isentropic  layers  and  vertical  redis¬ 
tribution  of  condensed-phase  HN03  and  H20. 

[13]  The  LaRC  LCTM  is  initialized  using  03,  H20,  NOx 
(NO  +  N02),  CH4,  HC1,  HF,  and  aerosol  surface  area 
density  (SAD)  observations  from  HALOE  and  03,  H20, 
and  SAD  observations  from  POAM  III.  Species  considered 
in  the  model  but  not  observed  from  the  satellites,  termed 
auxiliary  species,  are  initialized  using  the  IMPACT  model 
climatology.  Auxiliary  HALOE  species  are  initialized  using 
HALOE  CH4  and  sunrise  and  sunset  regressions  between 
the  auxiliary  species  and  CH4  obtained  from  the  IMPACT 
model  climatology.  Auxiliary  POAM  III  species  are 
obtained  using  UKMO  PV  (at  the  POAM  III  sounding) 
and  a  sunrise  or  sunset  regression  between  the  auxiliary 
species  and  PV  from  the  IMPACT  model  climatology.  Only 
those  HALOE  observations  north  of  30°S  with  less  than 
30%  measurement  uncertainty  for  03,  CH4,  HF,  H20,  HC1, 
and  Northern  Hemisphere  POAM  III  observations  with  less 
than  30%  measurement  uncertainty  for  03  and  H20  are 
included  in  the  simulations.  The  Cly  is  adjusted  when  the 
sum  of  the  observed  HC1  and  auxiliary  CLN03  +  HOCL  is 
greater  than  the  auxiliary  Cly  and  the  observed  NO  and  N02 
are  adjusted  so  that  the  sum  of  the  NOy  species  does  not 
exceed  the  auxiliary  NOy.  Initialization  concludes  with  a  5- 
day  diurnal  adjustment  in  which  HALOE  03,  H20,  and  HC1 
and  POAM  III  03  and  H20  are  reset  to  their  initial  values 
after  each  time  step.  This  final  initialization  step  reduces 
the  large  initial  radical  tendencies  associated  with  chem¬ 
ical  imbalances  between  the  observations  and  IMPACT 
climatology. 

[14]  A  Lagrangian  parameterization  of  diffusive  mixing  is 
also  included  in  the  LaRC  LCTM.  This  parameterization  is 
based  on  n-member  mixing  concepts  and  has  a  mixing 
efficiency  factor  that  is  a  function  of  the  relative  strength 
of  strain  and  rotation  in  the  flow  (Q)  [Fairlie  et  al,  1999]. 
Trajectories  within  regions  of  predominantly  rotational  flow 
have  negative  Q  and  mix  with  neighboring  trajectories  on  a 
weak  background  mixing  timescale  of  60  days.  Trajectories 
within  regions  of  shear  have  positive  Q  and  the  mixing 
timescales  can  be  as  fast  as  1-2  days  under  strong  shear 
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deformation.  The  flow  dependent  diffusion  results  in  epi¬ 
sodic  enhancements  in  mixing  efficiencies  for  a  particular 
trajectory  as  it  moves  in  and  out  of  shear  zones.  Sensitivity 
tests  have  shown  that  the  mixing  parameterization  performs 
well  in  both  strong  and  weak  mixing  regimes.  See  Fairlie  et 
al  [1999]  for  a  complete  discussion  of  the  physical  basis  for 
the  diffusion  parameterization. 

[15]  The  time  constants  for  mixing  are  the  same  as  given 
by  Fairlie  et  al  [1999].  However,  the  original  formulation 
was  found  to  be  globally  nonconservative.  The  nonconser¬ 
vation  arose  from  two  assumptions  that  changed  the  mean 
mixing  ratio  of  the  trajectories  (n-members)  involved  in 
each  mixing  event.  First,  mixing  tendencies  were  evaluated 
using  different  efficiencies  for  each  trajectory,  resulting  in 
more  rapid  approach  toward  the  n-member  mean  for  some 
trajectories  than  others  and  changes  in  the  n-member  mean. 
Second,  a  fixed  radius  was  used  to  determine  the  n-mem¬ 
bers  involved  in  mixing  with  each  trajectory.  Consequently, 
a  particular  trajectory  could  be  involved  in  more  than  1  n- 
member  mixing  event  (overlap),  resulting  in  changes  in  the 
combined  n-member  mean.  To  assure  global  conservation 
we  now  choose  the  n-member  mixing  ensemble  based  on 
the  UKMO  assimilation  grid  resolution  (all  trajectories 
within  1  vertical  grid  box,  2  north-south,  and  2  east-west 
neighboring  grid  boxes  mix),  thus  avoiding  overlap.  To 
assure  conservation  of  the  mean  mixing  ratio  during  each 
mixing  event  we  now  use  a  9-point  averaged  Q  from  the 
UKMO  grid  to  determine  the  mixing  efficiency  for  all 
trajectories  in  the  event.  The  new  formulation  results  in 
episodic  enhancements  in  the  mean  absolute  amplitude  of 
the  diffusive  tendencies  within  the  Arctic  vortex  as  internal 
shear  zones  develop  in  response  to  vortex  deformation,  but 
cannot  instantaneously  change  the  mean  mixing  ratio  within 
the  vortex  except  by  mixing  across  the  vortex  edge. 

[16]  The  vortex  edge  definition  used  in  the  model  simu¬ 
lations  is  closely  related  to  the  “kinematic  edge”  defined  by 
Pierce  and  Fairlie  [1993].  The  vortex  edge  is  determined  by 
integrating  Q  and  vorticity  around  rotational  streamlines  of 
the  flow.  The  northernmost  streamline  with  negative  inte¬ 
grated  Q  and  positive  (cyclonic)  integrated  vorticity  defines 
the  edge  of  polar  vortex.  This  kinematic  edge  typically 
resides  on  the  polarward  flank  of  the  polar  night  jet  [Pierce 
et  al ,  1994],  is  similar  to  dynamical  vortex  edge  definitions 
such  as  defined  by  Nash  et  al  [1996]  but  is  associated  with 
a  minimum  in  the  stretching  rates  of  material  lines  (mini¬ 
mum  Liapunov  exponent)  instead  of  sharp  gradients  in 
potential  vorticity.  To  assure  that  the  physical  characteristics 
of  the  kinematic  boundary  are  accurately  represented  in  the 
model  we  do  not  allow  enhanced  mixing  to  occur  across  the 
edge  of  the  vortex.  Consequently,  the  diffusion  during  a 
particular  model  time  step  does  not  significantly  affect  the 
mean  mixing  ratios  within  the  vortex.  However,  as  trajecto¬ 
ries  descend  within  the  Arctic  vortex  they  mix  with  different 
populations  at  different  altitudes.  Consequently,  the  cou¬ 
pling  between  diabatic  descent  and  mixing  can  change  the 
mean  mixing  ratio  within  the  vortex  at  a  given  altitude. 

3.  The  Solve  LaRC  LCTM  Simulation 

[17]  The  model  was  integrated  from  1  October  1999  to  28 
March  2000.  Trajectories  from  1995  HALOE  and  1312 
POAM  III  occultations  were  included  in  the  model  integra¬ 


tion.  Individual  trajectories  were  allowed  to  exist  for  60 
days.  This  long  lifetime  of  the  trajectories  introduces  errors 
in  the  trajectory  position  (particularly  longitude)  but  by 
considering  the  ensemble  mean  of  a  large  number  of 
trajectories  the  effects  of  these  errors  are  minimized  [Pierce 
et  al,  1999].  Figure  1  shows  the  sampling  characteristics  of 
the  HALOE  and  POAM  III  occultation  measurements  over 
the  course  of  the  1999/2000  winter  along  with  zonal  mean 
wind  speeds  at  450K.  In  the  October-November  time 
period  (Julian  days  —90  through  —30)  HALOE  sampled 
extensively  equatorward  of  60°N  with  measurements  within 
the  vortex  during  the  northern  most  extent  of  the  orbit. 
POAM  III  measurements  became  available  at  the  beginning 
of  November  (Julian  day  —60)  and  sampled  poleward  of 
60°N.  HALOE  experienced  data  losses  due  to  elevated 
instrument  temperatures  and  limited  spacecraft  power  dur¬ 
ing  much  of  November  1999.  The  Arctic  vortex  was  just 
beginning  to  form  at  this  time  with  wind  speeds  generally 
near  15  m/s  at  450K.  During  the  December- January  time 
period  (Julian  days  —30  through  30)  POAM  III  sampled 
generally  poleward  or  near  the  core  of  the  mean  polar  night 
jet  at  45 OK  while  HALOE  provided  primarily  midlatitude 
measurements  to  the  south  of  the  jet,  which  strengthened 
during  this  period,  reaching  mean  wind  speeds  up  to  32  m/s 
in  January.  HALOE  experienced  additional  data  losses 
during  late  January  and  early  February.  During  the  Febru¬ 
ary -March  time  period  (Julian  days  30  through  90)  the 
Arctic  vortex  showed  more  transient  behavior  with  a  rapid 
deceleration  of  the  zonal  mean  45 OK  jet  during  the  last  half 
of  February.  POAM  III  and  HALOE  sampled  both  inside 
and  outside  the  Arctic  vortex  during  this  time  period. 

[is]  During  the  SOLVE  deployments  we  found  that  the 
ensemble  trajectories  initialized  from  the  HALOE  and 
POAM  III  measurements  did  not  adequately  sample  the 
cold  pool  within  the  interior  of  the  Arctic  vortex  in  the 
lower  stratosphere  (ensemble  averaged  trajectory  temper¬ 
atures  were  considerably  warmer  than  the  mean  vortex 
temperatures  obtained  from  the  UKMO  assimilation).  To 
alleviate  this  sampling  bias  we  developed  a  cold  pool 
initialization  procedure,  which  involves  “spawning”  of 
new  “daughter”  trajectories  near  the  “parent”  trajectory. 
Each  daughter  has  the  same  chemical  characteristics  as  the 
parent  and  was  initialized  100  m  in  front  of  its  parent  along 
the  trajectory  path.  Shear  deformation  leads  to  the  dispersal 
of  the  daughter  trajectories  with  time.  Conditions  for 
spawning  were  evaluated  once  every  2  days  (at  12Z) 
throughout  the  integration  and  spawning  occurred  whenever 
the  parent  trajectory  age  was  greater  than  30  days,  the 
latitude  was  greater  than  30°N,  and  the  temperature  was 
less  than  192.5K.  As  will  be  shown,  the  cold  pool  initial¬ 
ization  resulted  in  improved  sampling  of  the  cool  pool.  The 
spawning  procedure  means  that  a  particular  parent  trajectory 
can  potentially  spawn  a  number  of  different  daughter 
trajectories,  particularly  if  it  remains  within  the  cold  pool 
for  a  number  of  days,  effectively  extending  its  lifetime. 
However,  no  new  chemical  information  is  being  added  to 
the  simulation  and  the  cold  pool  initialization  merely 
improves  the  sampling  characteristics  of  the  ensemble 
trajectory  simulation. 

[19]  The  first  60  days  of  the  simulation  (October- 
November  1999)  are  considered  a  “spin-up”  period  during 
which  the  Arctic  vortex  forms  and  the  HALOE  and 
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Figure  1.  Sampling  characteristics  of  the  HALOE  and  POAM  occultation  measurements  between 
October  1999  (Julian  day  —90)  and  March  2000  (Julian  day  90).  Circles  are  POAM  III  occultation  points 
and  triangles  are  HALOE  occultation  points.  The  zonal  mean  wind  speed  (m/s)  on  the  45 OK  potential 
temperature  surface  is  contoured. 


POAM  III  trajectories  accumulate.  The  rest  of  this  paper 
will  focus  on  the  December- March  1999/2000  portion  of 
the  simulation. 

3.1.  Comparisons  With  Early  Winter  Balloon  FTIR 
Measurements 

[20]  On  3  December  1999  the  balloonbome  MkIV  Four¬ 
ier  transform  infrared  (FTIR)  solar  absorption  spectrometer 
[Toon,  1991]  was  launched  into  the  core  of  the  Arctic  vortex 
(equivalent  latitude  of  80°N).  The  intent  of  this  flight  was  to 
characterize  the  state  of  the  Arctic  vortex  in  early  winter. 
Comparisons  between  the  single  MkIV  profile  and  the  range 
of  mixing  ratios  predicted  by  the  LaRC  LCTM  provides  an 
assessment  of  how  well  the  MkIV  profile  represents  the 
vortex  as  a  whole.  The  MkIV  measurements  also  provide  an 
opportunity  to  assess  the  LaRC  LCTM  predictions  during 
the  early  winter,  and  establish  how  well  the  model  has 
characterized  the  Arctic  vortex  prior  to  the  period  of 
enhanced  photochemical  loss.  Figure  2  shows  observed 
profiles  of  03,  NOy,  and  Cly  as  a  function  of  potential 
temperature,  along  with  the  individual  trajectory  predictions 
of  these  species  concentrations  within  the  vortex,  and  the 
predicted  vortex  ensemble  mean  profiles.  A  synoptic  map  of 
the  distribution  of  the  ensemble  ozone  prediction  at  575K- 
625K  is  also  shown.  The  bold  black  line  in  the  synoptic  map 
shows  the  edge  of  the  Arctic  vortex  at  600K,  the  thick,  short 


black  line  over  Sweden  is  the  trajectory  of  the  MkIV,  and 
the  contours  indicate  the  wind  speed  at  600K. 

[21]  The  observed  03  profile  falls  within  the  range  of 
ensemble  ozone  mixing  ratios  predicted  by  the  model  but  is 
clearly  on  the  lower  side  of  the  model  prediction,  which 
shows  a  wide  range  of  03  mixing  ratios  within  the  vortex 
above  approximately  500K.  Trajectories  initialized  with 
HALOE  and  POAM  III  observations  show  a  similar  spread 
in  03  mixing  ratios  although  those  initialized  with  HALOE 
data  tend  to  show  a  bimodal  distribution  above  800K  while 
those  initialized  with  POAM  III  observations  are  more 
uniformly  distributed  between  high  and  low  mixing  ratios. 
The  origin  of  the  bimodal  HALOE  03  is  not  immediately 
clear.  Investigation  of  the  initial  characteristics  of  these 
trajectories  establishes  that  it  is  not  associated  with  the 
latitude,  time,  or  altitude  at  which  the  HALOE  measurement 
occurred.  However,  the  synoptically  mapped  ensemble  03 
distribution  shows  that  the  low  03  mixing  ratios  are  found 
in  the  interior  of  the  Arctic  vortex,  well  poleward  of  the 
polar  night  jet,  while  the  higher  03  mixing  ratios  are  found 
near  the  edge  of  the  vortex  at  these  potential  temperatures. 

[22]  The  predicted  vortex  ensemble  mean  NOy  and  Cly 
profiles  are  very  close  the  observed  MkIV  profiles  between 
400K  and  900K  with  mean  absolute  differences  of  1 1  %  and 
6%,  respectively  indicating  that  mixing  ratios  for  the 
chlorine  and  nitrogen  reservoirs  are  well  characterized  by 
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Figure  2.  Vertical  profiles  of  03,  NOy,  and  Cly  observed  by  the  MkIV  FTIR  instrument  and  predicted 
by  the  model  on  3  December  1999  (12/03).  The  observed  profile  is  indicated  by  the  line  with  diamonds. 
The  predicted  vortex  ensemble  mean  profile  is  indicated  by  the  line  with  triangles.  Individual  vortex 
trajectory  mixing  ratios  are  indicated  by  red  circles  for  trajectories  initialized  from  POAM  III  occultations 
and  blue  triangles  for  trajectories  initialized  from  HALOE  occultations.  A  synoptic  map  of  the 
distribution  of  the  ensemble  ozone  prediction  (ppmv)  at  575K-625K  is  shown  in  the  lower  right  panel. 
The  bold  black  line  in  the  synoptic  map  shows  the  edge  of  the  Arctic  vortex  at  600K,  the  thick,  short 
black  line  over  Sweden  is  the  trajectory  of  the  MkIV.  The  contours  indicate  the  wind  speed  (m/s)  at  600K. 


the  PV  and  CH4  mapped  IMPACT  model  climatology  at 
these  altitudes.  The  comparison  between  the  LaRC  LCTM 
ensemble  NOy  and  Cly  predictions  and  MkIV  observations 
indicate  that  the  observed  profile  provides  a  reasonable 
characterization  of  the  mean  characteristics  of  the  early 
1999/2000  winter  vortex  below  800K. 

[23]  The  LaRC  LCTM  predictions  indicate  that  there  was 
a  radial  gradient  in  the  03  mixing  ratios  within  the  Arctic 
vortex  at  600K.  A  similar  distribution  of  03  within  the  early 
winter  Arctic  vortex  was  found  by  Kawa  et  al.  (unpublished 
manuscript,  2001).  However,  our  vortex  edge  definition 


suggests  that  the  03  gradient  begins  within  the  vortex  while 
Kawa  et  al.  found  that  the  gradient  lies  between  the  inner 
and  outer  vortex  edge  defined  by  Nash  et  al  [1996].  This 
apparent  discrepancy  arises  due  to  particularly  large  differ¬ 
ences  between  the  Nash  vortex  edge  and  the  kinematic  edge 
on  this  day.  On  average,  the  differences  between  the  kine¬ 
matic  edge  and  the  Nash  et  al  [1996]  edge  are  negligible. 
Comparisons  between  the  Nash  et  al  [1996]  and  kinematic 
edge  during  the  SOLVE  time  period  (not  shown)  indicate 
that  mean  differences  are  less  than  2°  in  equivalent  latitude, 
with  the  kinematic  edge  systematically  equatorward  of  the 
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Nash  edge.  Consequently,  the  kinematic  edge  is  roughly 
equivalent  to  the  “outer”  edge  as  defined  by  Nash  et  al. 
[1996].  What  is  clear  from  Figure  2  is  that  the  radial 
gradient  in  03  at  600K  is  poleward  of  the  polar  night  jet 
core.  As  will  be  shown  later,  air  from  500K  to  600K  is 
advected  poleward  and  then  descends  to  near  45 OK  by 
spring  where  it  experiences  enhanced  photochemical  loss. 
Because  of  this  descent,  the  03  mixing  ratios  between  500K 
and  600K  during  early  winter  play  a  critical  role  in 
observational  determinations  of  the  photochemical  loss 
rates.  The  range  of  03  mixing  ratios  predicted  within  the 
vortex  at  600K  during  the  beginning  of  December  compli¬ 
cates  the  diagnosis  of  springtime  photochemical  ozone  loss 
in  the  lower  stratosphere  from  observed  03-tracer  correla¬ 
tions.  The  LaRC  LCTM  predictions,  which  are  initialized 
with  observed  03,  suggest  that  assuming  that  a  few  early 
winter  03  profiles  can  characterize  the  “initial”  state  of  the 
Arctic  vortex  may  be  an  unrealistic  assumption. 

[24]  A  number  of  different  photochemical,  microphysical, 
and  diffusive  processes  must  be  simulated  along  each 
trajectory  to  accurately  represent  the  mean  vortex  evolution. 
Figure  3  illustrates  these  processes  for  a  typical  vortex 
trajectory.  This  particular  trajectory  was  spawned  from  a 
trajectory  originating  from  a  POAM  III  observation  at 
525K.  The  trajectory  was  initialized  in  the  beginning  of 
December  1999  at  517K  and  descended  to  406K  during  its 
lifetime.  It  remained  within  the  vortex  until  late  March, 
making  numerous  orbits  around  the  Arctic  vortex  and 
encountering  temperatures  below  195K  for  extended  peri¬ 
ods  of  time  between  late  December  and  mid-January  (Julian 
days  —10  to  20).  Two  periods  were  cold  enough  for  ice 
formation.  The  trajectory  spent  most  of  December  and 
January  orbiting  in  and  out  of  the  region  of  continuous 
polar  night  with  typical  daylight  exposure  times  of  6  hrs  per 
day  (0.25  daylight  fraction).  By  late  January  (Julian  day  30) 
the  daylight  exposure  of  the  trajectory  began  to  increase 
while  the  03  mixing  ratio  of  the  trajectory  begins  to  show 
more  rapid  declines.  Diffusive  03  tendencies  can  be  quite 
large  but  very  episodic  as  the  trajectory  moves  in  and  out  of 
regions  of  enhanced  mixing  within  the  vortex.  Diffusive 
tendencies  account  for  much  of  the  high  frequency  variation 
in  the  03  mixing  ratio  along  the  trajectory.  Enhanced 
photochemical  ozone  loss  occurs  during  periods  of  higher 
daylight  fraction.  During  December  through  mid-January, 
the  peak  net  loss  rates  (production-loss:  03  P-L  in  Figure  3) 
for  this  particular  trajectory  range  from  20  to  40  ppbv/day. 
By  late  January  the  net  loss  rates  accelerate  with  peak  loss 
rates  near  70  ppbv/day  occurring  at  the  beginning  of  March 
(Julian  day  60).  Chlorine  activation,  diagnosed  from  the 
Clx/Cly  ratio,  begins  in  early  December  and  reaches  nearly 
80%  by  mid- January.  This  particular  trajectory  shows 
chlorine  deactivation  from  late  January  through  February 
and  then  a  period  of  reactivation  in  March  associated  with 
colder  temperatures  near  Julian  day  60. 

[25]  The  trajectory  NOy  shows  a  series  of  small  declines 
during  December,  prior  to  the  first  encounter  with  temper¬ 
atures  below  ice  frost  points.  These  declines  are  associated 
with  the  formation  of  condensed-phase  HN03  and  subse¬ 
quent  sedimentation,  which  removes  HN03  from  the  tra¬ 
jectory  NOy  at  about  0.3  ppbv/day.  The  first  ice  event  for 
this  trajectory  occurs  during  the  last  week  of  December 
(Julian  day  —5)  with  declines  in  NOy  of  approximately  3 


ppbv/day,  associated  with  sedimentation  of  HN03.  This  ice 
event  occurs  during  a  period  of  high  diffusive  tendencies 
that  offset  much  of  the  denitrification  tendencies.  The 
largest  change  in  NOy  is  associated  with  a  strong  diffusive 
mixing  event  in  the  beginning  of  January  that  results  in 
significant  reductions  in  the  trajectory  NOy.  This  mixing 
event  is  followed  by  a  second  ice  event  that  occurs  in  early 
January  (Julian  day  10)  and  brings  the  NOy  mixing  ratio  of 
the  trajectory  to  nearly  zero.  This  complex  interaction 
between  sedimentation  and  mixing  occurs  as  diabatic 
descent  and  localized  sedimentation  leads  to  the  develop¬ 
ment  of  sharp  NOy  gradients  within  the  vortex  and  neigh¬ 
boring  trajectories  with  a  wide  range  of  NOy  mixing  ratios 
mix. 

[26]  It  is  clear  from  the  evolution  of  this  single  vortex 
trajectory  that  each  trajectory  encounters  a  broad  range  of 
chemical  and  dynamical  environments  over  the  course  of 
the  winter.  To  accurately  characterize  the  evolution  of  the 
Arctic  vortex  one  must  look  at  the  statistics  of  a  number  of 
trajectories  to  determine  the  ensemble  mean  state  of  the 
vortex.  Ensemble  averaging  is  accomplished  by  daily 
binning  all  vortex  trajectories  in  potential  temperature  bins 
defined  by  the  isentropic  database  and  averaging  to 
determine  daily  profiles  of  the  vortex  ensemble  mean 
chemical  and  dynamical  characteristics.  Time  series  of 
the  vortex  ensemble  mean  evolution  are  presented  in  the 
next  section. 

3.2.  Vortex  Ensemble  Mean  Evolution 

[27]  Figure  4  shows  the  time  series  of  the  vortex  aver¬ 
aged  temperatures  obtained  from  the  gridded  UKMO 
temperatures,  as  well  as  ensemble  averages  of  initial 
potential  temperature,  condensed  HN03  (NAT),  and  tem¬ 
perature  for  all  vortex  trajectories  as  a  function  of  potential 
temperature.  The  Arctic  cold  pool  forms  above  500K 
during  the  later  half  of  December  1999.  The  cold  pool 
gradually  descends  over  the  course  of  the  winter  with 
vortex  averaged  UKMO  temperatures  remaining  below 
202K  until  the  end  of  January  2000.  Manney  et  al. 
[2002]  compared  1999/2000  Arctic  temperatures  from  a 
number  of  different  assimilation  centers  including  UKMO. 
They  found  that  the  UKMO  cold  pool  was  significantly 
larger  than  found  other  assimilation  data  sets  during 
January  2000,  and  that  the  UKMO  coldest  temperatures 
where  generally  colder  than  other  assimilated  temperatures. 
However,  trajectory  studies  conducted  by  Manney  et  al 
[2002]  indicated  that  the  UKMO  temperature  cold  bias  did 
not  result  in  significant  bias  in  the  extent  of  PSC  exposure. 
During  the  beginning  of  February  there  is  a  minor  dis¬ 
turbance  that  results  in  a  temporary  warming  of  the  vortex, 
but  colder  temperatures  return  by  mid-February  and  persist 
until  mid-March,  when  the  final  warming  begins  and  the 
vortex  temperatures  warm  rapidly.  Diabatic  descent  is 
evident  in  the  ensemble  mean  initial  potential  temperature 
of  the  vortex  trajectories.  Trajectories  that  originated  at  550K 
are  found  at  475K  during  early  December  and  descend  to 
45 OK  by  late  February  and  early  March.  The  550K  isentrope 
is  near  the  lower  limit  of  altitudes  within  the  vortex  with  a 
wide  range  of  03  mixing  ratios  in  early  December  (see 
Figure  2).  The  overturning  of  the  initial  isentropes  during 
mid-February  is  a  consequence  of  intrusions  of  middle 
latitude  air  into  the  vortex  at  500K.  The  ensemble  mean 


SOL  60  -  8 


PIERCE  ET  AL.:  CHEMICAL  EVOLUTION  OF  THE  ARCTIC  VORTEX 


Temperature 


Time  (Julian  Day) 


Daylit  Fraction  03 


Time  (Julian  Day)  Time  (Julian  Day) 


03  Diffusion 


Time  (Julian  Day) 


MOy 


Time  (Julian  Day) 


03  P-L 


Time  (Julian  Day) 

MOy  Diffusion 


Time  (Julian  Day) 


Clx/Cly  Ratio 


-40  20  0  20  40  60  60 

Time  (Julian  Day) 


NOy  Evaporation+Sedimentation 


Time  (Julian  Day) 


Figure  3.  Time  series  of  photochemical,  microphysical,  and  diffusive  processes  along  a  typical  vortex 
trajectory.  Panels  in  the  upper  row  show  trajectory  temperature  (K),  daylight  fraction,  and  03  mixing  ratio 
(ppmv).  Panels  in  the  middle  row  show  diffusive  and  photochemical  tendencies  (ppbv/day),  and  the  ratio 
of  Clx/Cly  for  the  trajectory.  Panels  in  the  bottom  row  show  the  trajectory  NOy  mixing  ratio  (ppbv),  and 
diffusive  and  evaporation  +  sedimentation  tendencies  (ppbv/day).  The  vertical  dotted  lines  indicate  days 
which  ice  PSCs  were  encountered  by  the  trajectory.  Time  is  indicated  by  Julian  day  beginning  20 
November  (Julian  day  —40). 


condensed-phase  HN03  indicates  the  presence  of  PSCs 
within  the  vortex.  PSCs  initially  form  at  the  beginning  of 
December  between  500K  and  5  5 OK  and  grow  in-depth  to 
between  400K  and  600K  (15.5  and  23.5  km)  by  late 
December.  The  PSCs  persist  within  the  vortex  until  early 
February,  when  warming  of  the  vortex  results  in  evaporation 


of  the  NAT  particles.  Limited  PSCs  reappear  during  late 
February  and  early  March  between  400K  and  450K. 

[28]  The  ensemble  mean  temperatures  from  the  LaRC 
LCTM  trajectories  capture  the  vertical  and  temporal  varia¬ 
bility  found  in  the  vortex  averaged  UKMO  gridded  temper¬ 
atures  very  well.  This  agreement  between  the  evolution  of  the 
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Figure  4.  Potential  temperature-time  series  of  the  vortex  averaged  temperatures  (K)  obtained  from  the 
gridded  UKMO  temperatures  (upper  panel)  and  vortex  ensemble  means  of  initial  potential  temperature 
(ThO),  condensed  HN03  (NAT)  (ppbv),  and  temperature  (tmp)  (K).  The  first  day  of  each  month  is 
indicated. 


ensemble  mean  temperatures  and  the  evolution  of  the  vortex 
averaged  UKMO  gridded  temperatures  indicates  that  the 
LaRC  LCTM  is  adequately  sampling  the  meteorological 
conditions  within  the  Arctic  vortex  during  this  period.  How¬ 
ever,  the  coldest  ensemble  mean  vortex  temperatures  are  2K 
colder  then  the  coldest  vortex  averaged  UKMO  temperatures 
during  mid- January  and  late  February.  The  colder  ensemble 
mean  temperatures  during  this  period  arise  because  of  the 
dense  sampling  within  the  Arctic  cold  pool  associated  with 
the  cold  pool  initialization  procedure  discussed  previously. 
This  sampling,  coupled  with  the  general  cold  bias  in  the 
UKMO  temperatures  (S.  Buss  et  al.,  Arctic  stratospheric 
temperature  in  the  winters  1999/2000  and  2000/2001:  A 
quantitative  assessment  and  microphysical  implications, 


submitted  to  Journal  of  Geophysical  Research,  2001)  will 
tend  to  bias  the  LaRC  LCTM  results  toward  a  more  exten¬ 
sively  processed  vortex  than  may  have  actually  occurred. 

[29]  Figure  5  shows  time  series  of  ensemble  mean  vortex 
NOy,  and  NOy  tendencies  due  to  evaporation-sedimentation, 
diffusion,  and  vertical  advection.  The  ensemble  mean  abso¬ 
lute  amplitude  of  the  diffusive  tendencies  is  shown  to 
identify  time  periods  where  diffusive  tendencies  are  likely 
to  play  a  significant  role  in  determining  the  constituent 
distributions  within  the  vortex.  (Recall  that  diffusive  ten¬ 
dencies  have  negligible  effect  on  the  vortex  ensemble 
mean).  The  ensemble  mean  vertical  advection  tendency  is 
obtained  by  determining  the  ensemble  mean  diabatic 
descent  rate  and  multiplying  it  by  the  vertical  gradient  of 
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Figure  5.  Potential  temperature-time  series  of  vortex  ensemble  mean  NOy  (ppbv),  NOy  evaporation 
plus  sedimentation  tendency  (pptv/day),  absolute  value  of  NOy  diffusive  tendency  (pptv/day),  and  NOy 
vertical  advection  tendency  (pptv/day).  Negative  tendencies  are  indicated  by  dashed  contours.  The  first 
day  of  each  month  is  indicated. 


the  ensemble  mean  vortex  NOy  using  potential  temperature 
as  a  vertical  coordinate. 

[30]  NOy  isopleths  gradually  descend  due  to  vertical 
advection  during  the  first  half  of  December  1999  with 
negligible  contributions  from  evaporation,  sedimentation, 
or  diffusion.  By  late  December  the  extent  of  cold  temper¬ 
atures  becomes  large  enough  for  significant  sedimentation 
of  condensed  HN03  to  occur  and  the  vortex  is  denitrified  at 
a  rate  of  over  400  pptv/day.  This  rapid  denitrification  occurs 
in  conjunction  with  dehydration  that  reduces  the  ensemble 
mean  H20  by  30  ppbv/day  during  the  same  period  (not 
shown).  Comparison  with  Figure  4  (TH0  and  condensed 
HN03)  shows  that  this  period  of  intense  sedimentation 
corresponds  to  a  period  where  the  isopleths  of  condensed 
HN03  descend  more  rapidly  than  the  initial  potential 


temperature  isopleths  between  500K  and  45  OK,  indicating 
cross-isentropic  exchange  of  condensed-phase  HN03  via 
evaporation  of  falling  particles.  Vertical  advection  of  higher 
NOy  continuously  resupplies  NOy  at  a  rate  of  150  ppbv/day 
near  the  top  of  the  cold  pool.  Local  NOy  losses  due  to 
denitrification  overwhelm  the  ability  of  vertical  advection  to 
resupply  NOy  within  the  cold  pool  as  is  clearly  evident  in 
the  rapid  decline  in  ensemble  averaged  NOy  during  January. 
This  period  of  rapid  denitrification  terminates  by  mid- 
January  leaving  approximately  5-6  ppbv  of  NOy  at  45 OK 
throughout  the  rest  of  the  winter.  Vertical  advection  con¬ 
tinues  to  bring  down  higher  NOy  from  aloft  at  altitudes 
between  45 OK  and  5 5 OK  while  vertical  advection  reduces 
NOy  mixing  ratios  above  550K.  The  ensemble  mean  abso¬ 
lute  amplitude  of  the  NOy  diffusive  tendency  shows  strong 
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Figure  6.  Potential  temperature-time  series  of  vortex  ensemble  mean  Cly  (ppbv),  Clx/Cly  ratio,  absolute 
value  of  Cly  diffusive  tendency  (pptv/day),  and  Cly  vertical  advection  tendency  (pptv/day).  Negative 
tendencies  are  indicated  by  dashed  contours.  The  first  day  of  each  month  in  indicated. 


mixing  within  the  vortex  during  the  period  of  rapid  deni¬ 
trification,  as  local  denitrification  and  vertical  advection 
leads  to  large  horizontal  gradients  in  NOy,  and  then  a  series 
of  diffusive  mixing  events  in  February  and  early  March. 
There  are  episodic  periods  where  the  ensemble  mean 
absolute  amplitude  of  the  diffusive  tendencies  are  compa¬ 
rable  with  the  tendencies  due  to  sedimentation. 

[31]  Figure  6  shows  time  series  of  ensemble  mean  vortex 
Cly,  Clx/Cly,  absolute  amplitude  of  diffusive  tendencies,  and 
vertical  advection  tendencies.  Throughout  most  of  the 
winter  the  Cly  isopleths  descend  due  to  vertical  advection. 
Cly  has  no  significant  loss  processes  in  the  lower  strato¬ 
sphere  but  the  partitioning  of  the  chlorine  species  within  Cly 
plays  a  critical  role  in  the  enhanced  ozone  loss  in  the  Arctic 
spring  [, Solomon ,  1999].  Chlorine  reservoir  species  (HC1 


and  CINO3)  are  converted  into  chlorine  radicals  (Clx  =  CIO 
+  2*C1202)  during  late  December  1999  and  throughout 
January  2000,  leading  to  ensemble  mean  Clx/Cly  ratios  in 
excess  of  0.8  between  450K  and  500K  in  late  January. 
Episodes  of  significant  diffusive  mixing  of  Cly  within  the 
vortex  occur  at  the  same  time  as  the  enhanced  NOy  mixing. 
Vertical  advection  increases  Cly  mixing  ratios  by  30  pptv/ 
day  between  500  and  45 OK  leading  to  a  very  sharp  vertical 
gradient  in  Cly  near  425K  by  February.  During  the  second 
half  of  February  negative  Cly  tendencies  due  to  vertical 
advection  dominate  above  approximately  45  OK,  leading  to 
reductions  in  the  ensemble  mean  vortex  Cly  mixing  ratios. 

[32]  Figure  7  shows  time  series  of  ensemble  mean  vortex 
03,  03  net  photochemical  loss  (P-L),  absolute  amplitudes  of 
diffusive  tendencies,  and  vertical  advection  tendencies. 
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Figure  7.  Potential  temperature-time  series  of  vortex  ensemble  mean  03  (ppmv),  03  photochemical 
production  -  loss  (P-L)  ( ppbv/day),  absolute  value  of  03  diffusive  tendency  ( ppbv/day),  and  03  vertical 
advection  tendency  (ppbv/day).  Negative  tendencies  are  indicated  by  dashed  contours.  The  first  day  of 
each  month  in  indicated. 


Vertical  advection,  with  tendencies  that  reach  30  ppbv/day 
at  450K,  largely  determines  the  evolution  of  03  in  the  lower 
stratospheric  Arctic  vortex  during  December  to  mid-January 
with  only  minor  contributions  from  photochemical  losses, 
which  are  near  5  ppbv/day  during  early  January.  There  is  a 
period  of  enhanced  photochemical  loss  between  450K  and 
500K  in  mid-February  but  the  major  loss  occurs  at  425K 
during  mid-March,  where  ensemble  mean  photochemical 
loss  rates  reach  over  60  ppbv/day.  These  rates  are  very  close 
to  the  peak  loss  rates  of  61.6  ±  4.8  ppbv/day  calculated  in 
mid-March  during  the  1999/2000  MATCH  campaign  [Rex 
et  al .,  2002].  The  ensemble  mean  absolute  amplitude  of  the 
03  diffusive  tendencies  shows  the  same  general  pattern  of 
episodes  of  strong  diffusion  above  45 OK  during  February 
and  March  as  found  in  NOy  and  Cly.  This  similarity  arises 


due  to  positive  vertical  gradients  in  all  three  constituents 
below  approximately  800K.  The  link  between  vertical 
constituent  gradients  and  diffusive  mixing  tendencies  will 
be  discussed  in  more  detail  later. 

[33]  Figure  8  shows  time  series  of  ensemble  mean  accu¬ 
mulated  photochemical  03  loss  and  the  relative  contribu¬ 
tions  from  the  C10x,  BrOx,  and  C10x  +  BrOx  loss  cycles  to 
the  total  photochemical  loss.  The  C10x  and  BrOx  loss  cycles 
are  defined  as  given  by  Wennberg  et  al.  [1994].  The 
accumulated  ensemble  mean  03  loss  by  the  end  of  March 
is  2.2  ppmv  at  45 OK.  This  accumulated  03  loss  is  slightly 
less  than  the  peak  accumulated  03  loss  of  2.7  ±  0.24  ppmv  at 
453K  obtained  by  the  MATCH  campaign  [Rex  et  al ,  2002]. 
Focusing  on  the  periods  of  significant  net  loss  reveals  that 
C10x  and  BrOx  loss  cycles  each  contribute  ^35%  to  the  net 
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Figure  8.  Potential  temperature-time  series  of  vortex  ensemble  mean  accumulative  03  photochemical 
production  -  loss  (P-L)  (ppmv),  fraction  of  loss  due  to  C10x  cycle,  fraction  of  loss  due  to  BrOx  cycle,  and 
fraction  of  loss  due  to  C10x  +  BrOx  cycles.  Net  loss  is  indicated  by  dashed  contours.  The  first  day  of  each 
month  in  indicated. 


loss  rates,  accounting  for  70%  of  the  net  loss  rates  during 
mid-February.  The  role  of  C10x  and  BrOx  loss  cycles 
increases  to  90-95%  of  the  net  loss  by  March  with  C10x 
cycles  accounting  for  55%  and  BrOx  cycles  accounting  for 
40-45%.  As  the  net  loss  rates  decline  in  mid  to  late  March 
the  relative  importance  of  BrOx  loss  cycles  increases  while 
the  importance  of  C10x  loss  declines.  This  shift  occurs  as 
chlorine  recovery  takes  place  and  CIO  mixing  ratios  decline. 
The  vortex  ensemble  mean  NOx  is  below  0.01  ppbv  (not 
shown)  during  early  March.  Which  compares  well  with  the 
mean  NOx  observed  by  instruments  onboard  the  ER2  within 
the  vortex  during  this  period  (0.02  ppbv).  Consequently,  the 
vortex  was  nearly  completely  denitrified  in  early  March 
2000  and  the  model  predicts  that  chlorine  recovery  occurs 
by  transport  of  extravortex  air  with  lower  CIO  mixing  ratios 
and  higher  N02  mixing  ratios  into  the  vortex. 


[34]  The  ensemble  mean  time  series  provide  a  useful 
measure  of  the  mean  characteristics  of  the  Arctic  vortex 
over  the  course  of  the  1999/2000  winter.  However,  they  do 
not  indicate  the  degree  of  variability  found  with  in  the 
vortex.  The  next  section  explores  the  vortex  variability 
during  early  March  (the  period  of  peak  photochemical 
losses)  in  detail  by  comparing  predicted  frequency  distri¬ 
butions  with  in  situ  observations  from  the  ER2.  These 
comparisons  also  serve  to  establish  the  validity  of  the 
ensemble  mean  predictions. 

3.3.  Vortex  Frequency  Distributions 

[35]  Figure  9  shows  frequency  distributions  (FDs)  of 
CH4,  03,  percent  denitrification,  and  Clx/Cly,  for  vortex 
measurements  and  coincident  model  trajectories  between 
425K  and  475K  taken  on  the  05/03,  07/03,  11/03,  and  12/03 
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Figure  9.  Frequency  distributions  of  CH4,  03,  percent  denitrification,  and  Clx/Cly,  for  vortex  in  situ 
measurements  (dotted  line)  and  coincident  model  trajectories  between  425K  and  475K  taken  on  the  05/ 
03,  07/03,  11/03,  and  12/03  ER2  flights.  Model  frequency  distributions  are  shown  for  simulations  with 
(solid  line)  and  without  (dashed  line)  diffusive  tendencies.  The  number  in  parenthesis  below  the  line 
caption  indicates  the  number  of  observations/trajectories  used  to  construct  the  frequency  distributions. 
The  coincidence  criteria  were  the  same  day,  10°  longitude,  5°  latitude,  and  5K  in  potential  temperature. 


ER2  flights.  The  coincidence  criteria  was  the  same  day,  10° 
longitude,  5°  latitude,  and  5K  in  potential  temperature, 
resulting  in  234  coincident  trajectories  for  these  flights. 
Model  results  include  a  simulation  where  diffusive  tenden¬ 
cies  were  turned  off  to  assess  the  role  of  diffusion  in 
determining  the  constituent  frequency  distributions.  The 


classification  of  ER2  measurements  as  vortex  measurements 
was  accomplished  by  interpolating  the  vortex  marker  field 
in  the  UKMO  isentropic  database  to  the  ER2  flight  track. 
3.3.1.  Ch4  Frequency  Distribution 
[36]  The  CH4  FD  provides  the  clearest  assessment  of  the 
impact  of  diffusion  on  the  simulation  since  it  is  long-lived  in 
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the  lower  stratosphere  and  its  vortex  FD  is  largely  deter¬ 
mined  by  the  interaction  between  vertical  descent  and  small- 
scale  diffusion.  The  predicted  vortex  CH4  FD  shows  a  strong 
peak  near  0.7  ppmv.  Which  is  in  closer  agreement  with  the 
observed  [Webster  et  al,  1994]  CH4  FD  (peak  near  0.9 
ppmv)  than  the  simulation  without  diffusive  mixing  which 
peaks  near  0.4  ppmv  and  has  a  broad  shoulder  extending  to 
1.0  ppmv.  The  differences  in  the  simulations  with  and 
without  diffusion  arise  due  to  differential  diabatic  descent 
(which  introduces  different  CH4  mixing  ratios  at  a  given 
altitude  due  to  the  vertical  gradient  in  CH4),  combined  with 
diffusive  mixing  which  tends  to  increase  the  low  CH4 
mixing  ratios  and  decrease  the  high  CH4  mixing  ratios.  This 
leads  to  more  uniform  CH4  mixing  ratios  (narrow  peak  in  the 
FD)  in  the  simulation  with  diffusion.  However,  the  model 
still  shows  a  bias  toward  low  CH4  mixing  ratios.  The  0.7 
ppmv  CH4  mixing  ratio  peak  found  in  the  predicted  FD  is 
close  to  the  ensemble  mean  vortex  CH4  mixing  ratio  (0.72 
ppmv)  at  570K  in  early  December  (not  shown),  which 
descends  at  a  rate  of  1.5K/day  during  December  through 
January  and  0.95K/day  during  February  reaching  450K  by 
early  March.  Greenblatt  et  al  [2002]  used  in  situ  and  remote 
balloon  measurements  to  assess  the  observed  descent  rates  of 
selected  CH4  isopleths  during  the  1999/2000  winter  and 
found  that  the  observed  0.72  ppmv  CH4  isopleth  descended 
from  5 5 OK  in  early  December  to  470K  in  early  March  at  a 
rates  of  l.OK/day  during  December  through  January  and 
0.6K/day  during  February.  The  larger  descent  rates  in  the 
model  CH4  isopleth  implies  stronger  diabatic  cooling  in  the 
vortex  than  actually  occurred.  In  addition  to  the  bias  intro¬ 
duced  by  stronger  diabatic  descent  there  is  a  sampling  bias 
introduced  by  the  cold  pool  spawning  procedure.  The 
trajectories  with  the  lowest  CH4  mixing  ratios  (near  0.4  in 
the  no-diffusion  FD)  encountered  the  cold  pool  in  late 
December  and  early  January  and  consequently  spawned 
daughter  trajectories.  The  spawning  leads  to  increased  num¬ 
bers  of  low  CH4  trajectories  that  introduce  a  bias  the  model 
FDs  toward  lower  CH4  mixing  ratios. 

[37]  The  modeled  CH4  FDs  have  a  few  CH4  mixing  ratios 
above  1.5.  These  high  CH4  mixing  ratios  suggest  that  there 
is  a  small  population  of  extravortex  trajectories.  These 
trajectories  must  have  entered  the  vortex  during  one  of  the 
periods  where  the  kinematic  barrier  was  weakened  by 
transient  disturbances,  allowing  horizontal  transport  into 
the  vortex  to  occur.  However,  the  trajectories  with  high 
CH4  mixing  ratios  make  up  a  very  small  fraction  of  the  total 
vortex  ensemble,  which  indicates  that,  based  on  the  model 
simulations,  horizontal  transport  of  lower  stratospheric, 
extravortex  air  into  the  vortex  was  relatively  weak  during 
the  1999/2000  winter. 

3.3.2.  03  Frequency  Distribution 

[38]  The  03  FDs  show  that  the  model  slightly  over¬ 
estimates  the  03  mixing  ratio  (peak  near  1.6  ppmv)  relative 
to  the  in  situ  measurements  [Proffitt  and  McLaughlin , 
1983],  which  show  a  peak  in  the  FD  near  1.5  ppmv  at 
these  altitudes,  and  has  a  broader  tail  at  high  03  mixing 
ratios  than  is  observed.  The  simulation  without  diffusion 
predicts  the  peak  at  1 .3  ppmv  and  the  distribution  is  clearly 
too  broad,  particularly  on  the  high  03  mixing  ratio  side  of 
the  peak.  The  tail  of  high  mixing  ratios  in  the  predicted  03 
could  be  due  to  over-sampling  of  higher  initial  mixing  ratios 
by  HALOE,  which  does  not  sample  the  deep  core  of  the 


Arctic  vortex  at  5  5 OK  during  the  early  winter  as  frequently 
as  POAM  III  does.  This  initial  high  bias  in  the  550K 
ensemble  prediction  is  evident  in  Figure  2,  which  shows 
that  trajectories  initialized  with  HALOE  03  mixing  ratios 
tend  to  cluster  near  the  high  03  portion  of  the  ensemble 
prediction  between  500K  and  800K. 

3.3.3.  Percent  Denitrification  Frequency  Distribution 

[39]  The  observed  and  predicted  percent  of  denitrification 
within  the  vortex  is  determined  by  subtracting  the  March  in 
situ  NOy  measurements  [Fahey  et  al ,  1989]  and  ensemble 
NOy  predictions  from  the  NOy  observed  during  the  12/03 
MkIV  balloon  flight,  using  N20  mixing  ratios  to  identify 
similar  air  masses.  The  03/05  in  situ  NOy  measurements  are 
filtered  to  exclude  PSC  encounters  since  condensed-phase 
HNO3  tends  to  introduce  uncertainties  in  the  gas-phase  NOy 
measurement.  No  PSCs  were  encountered  during  the  rest  of 
the  March  flights  considered  here.  The  predicted  and 
observed  distribution  of  denitrification  shows  relatively 
good  agreement  when  diffusion  is  included  although  the 
predicted  FD  is  somewhat  narrower  than  observed  with  a 
peak  near  60%  while  the  observed  peak  is  near  70% 
denitrification.  The  simulation  without  diffusion  shows  a 
broader  distribution  than  is  observed  with  a  peak  near  50%. 
The  observed  denitrification  shows  a  long  tail  extending 
toward  (and  beyond)  zero  that  is  not  present  in  the  predicted 
FD  and  is  too  broad  in  the  simulation  without  diffusion. 

[40]  The  10%  differences  between  the  location  of  the 
observed  and  predicted  peaks  in  the  percent  denitrification 
FDs  indicates  that  our  sedimentation  rates  may  be  slightly 
lower  than  actually  occur.  The  tail  in  the  observed  FD  could 
arise  due  to  re-evaporation  of  sedimenting  NAT  particles, 
incomplete  denitrification,  or  mixing  of  denitrified  air  with 
unprocessed  air  from  aloft.  Our  simulations  indicate  that  re¬ 
evaporation  of  sedimenting  NAT  contributes  approximately 
2  ppbv  of  NOy  at  these  altitudes  (shown  later).  Signatures  of 
mixing  with  unprocessed  air  from  aloft  are  evident  in  the 
predicted  NOy/CH4  correlations  (not  shown)  but  these  tra¬ 
jectories  tend  to  be  found  at  slightly  higher  potential  temper¬ 
atures  than  the  ER2  sampled  and  only  make  up  a  small 
fraction  of  the  vortex  ensemble  below  475K  during  March. 
Consequently,  it  is  most  likely  that  the  observed  tail  is 
associated  with  incomplete  denitrification  in  the  air  sampled 
by  the  ER2.  The  fact  that  the  no-diffusion  simulation  over¬ 
estimates  this  feature  while  the  simulation  with  diffusion 
underestimates  it  suggests  that  the  model  diffusion  may  be 
somewhat  stronger  than  actual  small-scale  mixing  processes. 

[41]  The  simple  bulk  equilibrium  parameterization  of 
denitrification  used  in  the  model  is  in  lieu  of  full  calcula¬ 
tions  involving  detailed  representation  of  the  particle  micro¬ 
physics  [Drdla  et  al ,  2002].  The  agreement  of  observed  and 
calculated  denitrification  indicates  that  the  model  parame¬ 
terization  captures  some  important  large-scale  features  of 
the  distribution  of  denitrified  air  in  the  Arctic  vortex. 
However,  observations  of  NOy,  HN03,  and  aerosol  particles 
during  SOLVE  have  revealed  the  presence  of  large  HN03- 
containing  particles  over  large  areas  of  the  vortex  in  the 
lower  stratosphere  [Fahey  et  al ,  2001;  Carslaw  et  al ,  2002; 
Northway  et  al ,  2002]  which  are  not  accounted  for  in  the 
model.  Particles  with  sizes  ranging  up  to  20  microns  in 
diameter  were  observed  throughout  late  January  to  mid- 
March  [Fahey  et  al ,  2001].  The  small  population  of  large 
particles  suggests  a  highly  selective,  nonequilibrium  nucle- 
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ation  process  must  be  responsible  for  the  formation  of  these 
particles  [Fahey  et  al .,  2001].  Instantaneous  flux  values 
indicate  that  these  particles  could  substantially  denitrify  the 
lower  stratosphere  in  a  matter  of  days.  In  contrast,  rapid 
denitrification  in  the  model  occurs  via  synoptic-scale  HN03 
removal  via  equilibrium  growth  of  a  NAT  (HNO3  3H20) 
coating  on  ice  particles. 

[42]  Davies  et  al  [2002]  conducted  chemical  transport 
model  simulations  during  the  1999/2000  winter  using  the 
SLIMCAT  model  to  explore  model  sensitivity  to  differences 
in  denitrification  parameterizations.  They  considered  two 
different  formulations,  an  ice  formation  scheme  very  similar 
to  our  scheme  and  a  NAT  scheme  that  assumed  a  bimodal 
size  distribution  (0.5  micron  and  6.5  micron  radii)  con¬ 
strained  by  the  Fahey  et  al  [2001]  observations.  In  the 
bimodal  NAT  scheme  the  small  mode  affects  heterogeneous 
rate  calculations  due  to  large  surface  area  densities  while  the 
large  mode  is  responsible  for  the  majority  of  the  denitrifi¬ 
cation  due  to  larger  fall  velocities.  Davies  et  al  [2002] 
found  that  both  schemes  produced  extensive  denitrification 
during  the  1999/2000  winter  when  UKMO  temperatures 
were  used  although  the  bimodal  NAT  scheme  produced 
somewhat  better  agreement  with  observed  denitrification 
than  the  ice  scheme.  In  particular,  the  ice  scheme  did  not 
reproduce  the  observed  extensive  denitrification  when 
ECMWF  temperatures  were  used. 

[43]  To  better  understand  the  role  of  ice  formation  in  the 
denitrification  processes  one  must  look  at  in  situ  water 
vapor  observations.  These  suggest  that  at  least  some  of 
the  observed  denitrification  could  have  occurred  in  associ¬ 
ation  with  dehydration.  In  situ  measurements  of  total  hydro¬ 
gen  at  the  highest  altitudes  sampled  by  the  ER2  (21km) 
suggest  that  widespread  dehydration  of  at  least  0.2  ppmv 
occurred  for  air  masses  with  CH4  mixing  ratios  of  0.7  ppmv, 
[Herman  et  al ,  2002],  and  balloon  borne  instruments 
indicate  dehydration  of  near  0.5  ppmv  at  CH4  mixing  ratios 
of  0.7  ppmv,  [Schiller  et  al ,  2002].  However,  in  situ 
observations  indicate  that  while  all  of  the  dehydrated  air 
was  denitrified,  much  of  the  denitrified  air  was  not  dehy¬ 
drated  (R.  Herman,  personal  communication,  2001).  The 
majority  of  the  simulated  denitrification  occurred  in  con¬ 
junction  with  dehydration. 

3.3.4.  Clx/Cly  Frequency  Distribution 

[44]  The  observed  FD  of  Clx/Cly  is  obtained  by  combin¬ 
ing  Cly  measurements  from  the  ACATS  IV  instrument 
[Elkins  et  al ,  1996],  with  CIO  +  2*CL202  measurements 
from  the  Harvard  C10x  instrument  (R.  Stimpfle  et  al., 
unpublished  manuscript,  2001).  The  observed  Clx/Cly  FD 
peaks  near  0.4  while  LaRC  LCTM  simulations  with  and 
without  diffusion  both  show  Clx/Cly  FDs  which  peak  near 
0.7 -0.9  with  long  tails  toward  low  Clx/Cly  ratios.  Diffusive 
mixing  doesn’t  significantly  change  the  partitioning  of  Clx 
within  the  Cly  family.  The  high  ratio  of  predicted  Clx/Cly 
relative  to  the  observed  ratio  is  larger  than  can  be  explained 
by  uncertainties  in  the  observed  Clx  (approximately  30%), 
and  indicates  that  the  model  chlorine  is  significantly  over 
activated  relative  to  the  measurements  during  early  March. 

3.3.5.  Chlorine  Frequency  Distributions 

[45]  To  further  explore  the  large  differences  between  the 
modeled  and  measured  chlorine  activation  levels  we  con¬ 
sider  the  partitioning  within  the  Cly  family.  Figure  10  shows 
FDs  of  observed  and  modeled  Clx,  Cly,  C1N03,  and  HC1. 


The  Cly  FDs  also  show  the  sum  of  the  observed  C1N03  + 
CIO  +  2C1202  +  HC1.  The  Clx  and  Cly  FDs  show  that  much 
of  the  bias  in  the  model  Clx/Cly  occurs  because  of  high 
model  Clx  since  the  modeled  Cly  is  in  excellent  agreement 
with  the  ACATS  measurements  when  diffusion  is  included. 
The  peak  in  the  sum  of  the  observed  chlorine  species  FD  is 
1.0  ppbv  lower  than  the  observed  Cly.  This  difference 
reflects  an  imbalance  in  the  observed  chlorine  budget  of 
approximately  33%.  The  response  of  the  model  Cly  FD  to 
diffusive  mixing  mirrors  the  CH4  FD  response  since  both 
are  conserved  in  the  lower  stratosphere  and  Cly  has  a 
positive  vertical  gradient  while  CH4  has  a  negative  vertical 
gradient.  The  excellent  agreement  between  the  modeled  and 
measured  Cly  occurs  in  spite  of  the  bias  in  the  modeled  CH4 
since  Cly  is  nearly  uniform  above  500K  (Figure  2). 

[46]  The  observed  FDs  of  the  chlorine  reservoir  species 
(CINO3  and  HC1)  show  peaks  near  0.4  and  0.3  ppbv 
respectively,  indicating  that  some  chlorine  recovery  has 
begun  for  the  majority  of  the  air  masses  observed  by  the 
ER2  during  early  March.  In  contrast,  the  model  C1N03  and 
HC1  FDs  show  peaks  near  zero,  indicating  that  most  of  the 
coincident  trajectories  are  still  highly  activated.  The  imbal¬ 
ance  in  the  observed  chlorine  budget  does  not  allow  a  clear 
assessment  of  the  model  Clx/Cly  bias  since  underestimates 
in  any  of  the  observed  chlorine  species  or  overestimates  in 
the  total  inorganic  chlorine  could  account  for  the  observed 
deficit.  Increases  the  observed  Clx  would  result  in  better 
agreement  between  the  modeled  and  measured  Clx/Cly, 
while  increases  in  the  observed  C1N03  or  HC1  would  result 
in  worse  agreement.  However,  the  FD  comparisons 
between  the  modeled  and  measured  chlorine  species  all 
indicate  that  the  model  chlorine  activation  is  high  relative  to 
the  measurements. 

4.  Synoptic  Case  Studies 

[47]  Results  from  the  ensemble  mean  time  series  and  FD 
analysis  raise  two  unanswered  questions.  The  first  question 
is  how  does  small-scale  diffusive  mixing  contribute  to 
changes  in  distribution  of  trace  gas  concentrations  within 
the  vortex?  The  second  question  is  why  does  the  model 
chlorine  remain  activated  during  early  March?  This  second 
question  is  particularly  critical  to  address  since  the  success 
in  predicting  the  observed  ozone  distribution  is  largely  do  to 
a  relatively  short  period  of  rapid,  chlorine  catalyzed  ozone 
losses  in  early  March,  just  prior  to  chlorine  deactivation, 
during  a  period  of  rapidly  changing  chlorine  partitioning 
and  significant  overestimates  in  modeled  Clx  relative  to  the 
measurements.  To  answer  these  questions  and  illustrate  the 
degree  of  spatial  variability  within  the  vortex  we  consider 
two  synoptic  case  studies. 

4.1.  Case  Study:  02/26/00 

[48]  The  first  case  study  considers  the  synoptic  situation 
during  the  ER2  flight  of  02/26/00.  This  flight  occurred 
during  one  of  the  episodes  of  large  ensemble  mean  ampli¬ 
tudes  in  diffusive  tendencies  (see  Figures  5,  6,  and  7). 
Figure  11  shows  synoptically  mapped  ensemble  predictions 
of  diabatic  heating  rates,  diffusive  mixing  efficiency  (Q), 
and  03  loss  rates  between  425K  and  475K.  A  comparison 
between  the  observed  and  predicted  correlation  between 
CH4  and  03  within  the  vortex  is  also  shown.  The  black  line 
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Coincident  Clx  Freq.  Dist.  425K-475K 


Coincident  Cly  Freq.  Dist.  425K-475K 


Figure  10.  Frequency  distributions  of  Clx,  Cly,  C1N03,  and  HC1,  for  vortex  in  situ  measurements 
(dotted  line)  and  coincident  model  trajectories  between  425K  and  475K  taken  on  the  03/05,  03/07, 
03/11,  and  03/12  ER2  flights.  Model  frequency  distributions  are  shown  for  simulations  with  (solid 
line)  and  without  (dashed  line)  diffusive  tendencies.  The  number  in  parenthesis  below  the  line  caption 
indicates  the  number  of  observations/trajectories  used  to  construct  the  frequency  distributions.  The 
coincidence  criteria  were  the  same  day,  10°  longitude,  5°  latitude,  and  5K  in  potential  temperature. 
The  Cly  frequency  distribution  (upper  right  panel)  also  shows  the  sum  of  the  observed  CINO3  + 
CIO  +  2C1202  +  HC1  (dash-dot  line). 
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425 K-  475K  Diabatic  Heating  Rate  02/26/00 


425K-  475K  Mixing  Efficiency  02/26/00 
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Figure  11.  Synoptically  mapped  ensemble  predictions  of  diabatic  heating  rates  (K/day),  diffusive 
mixing  efficiency  (Q)  (sec-1),  and  03  loss  rates  between  425K  and  475K  on  26  February  2000  (02/26/ 
00).  The  lower  right  panel  shows  a  comparison  between  the  observed  (black  triangles)  and  predicted 
(colored  crosses)  correlation  between  CH4  and  03  within  the  vortex.  The  black  line  in  the  synoptic  maps 
indicates  the  kinematic  edge  of  the  Arctic  vortex  at  45  OK.  The  bold  black  line  indicates  the  ER2  flight 
track. 


in  the  synoptic  maps  indicates  the  kinematic  edge  of  the 
Arctic  vortex  at  450K.  The  bold  black  line  indicates  the 
ER2  flight  track. 

[49]  Diabatic  heating  (time  rate  of  change  in  potential 
temperature  following  an  air  parcel)  is  the  vertical  velocity 
component  in  isentropic  coordinates  and  negative  diabatic 
heating  implies  descent  of  air  through  isentropic  surfaces. 
The  diabatic  heating  rates  for  the  ensemble  trajectories  show 
considerable  variability  within  the  vortex  with  maximum 
diabatic  cooling  (descent)  of  approximately  —  2.0K/day 
occurring  over  the  pole.  The  mixing  efficiency  shows  even 
more  variability  within  the  vortex  with  three  distinct  lobes 


of  strongly  negative  Q  (indicating  a  region  of  predominately 
rotational  flow  and  no  enhanced  diffusive  mixing)  and  a 
region  between  the  lobes  with  strongly  positive  Q  (indicat¬ 
ing  strong  shear  deformation  and  enhanced  diffusive  mix¬ 
ing).  The  shear  zones  and  rotational  lobes  develop  due  to 
distortions  in  the  Arctic  vortex  caused  by  midlatitude  anti- 
cyclonic  circulations  to  the  north  of  Alaska  and  central 
Siberia.  These  midlatitude  circulations  are  evident  in  the 
concave  curvature  of  the  kinematic  edge  and  weakly  neg¬ 
ative  ensemble  Q  in  these  regions. 

[50]  The  region  of  maximum  diabatic  descent  coincides 
with  this  enhanced  mixing  zone.  This  occurs  since  the 
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midlatitude  disturbances  which  cause  the  vortex  distortions 
also  transport  warm  air  poleward,  leading  to  enhanced 
radiative  cooling  near  the  polar  night  jet  core  [ Pierce  et 
al .,  1993].  The  largest  enhancements  in  diffusive  mixing 
tendencies  (not  shown)  occur  in  the  region  of  large  Q,  where 
both  positive  and  negative  tendencies  are  present.  The 
episodes  of  enhanced  diffusive  mixing  evident  in  Figures 
5,  6,  and  7  during  February  and  March  above  45 OK  arise 
because  of  the  coupling  between  the  vertical  constituent 
gradients,  differential  diabatic  descent,  and  small-scale  mix¬ 
ing  processes.  Strong  diabatic  descent  over  the  pole  intro¬ 
duces  a  population  of  trajectories  that  have  descended  from 
higher  altitudes  and  have  higher  03,  NOy,  and  Cly  mixing 
ratios  (and  lower  CH4  mixing  ratios)  than  their  neighboring 
trajectories,  into  the  region  of  enhanced  mixing.  This  pop¬ 
ulation  diffusively  mixes  with  the  neighboring  trajectories 
with  lower  (higher)  mixing  ratios.  Because  the  mean  mixing 
ratio  is  conserved  during  each  mixing  event,  the  mean 
mixing  ratios  of  the  two  populations  are  not  affected. 
However,  since  trajectories  encounter  different  populations 
at  different  altitudes  as  they  descend  within  the  vortex,  the 
mean  mixing  ratio  at  a  given  level  can  be  affected  by 
diffusion.  This  process  accounts  for  the  shifts  in  the  modeled 
frequency  distributions  when  diffusive  mixing  is  included  in 
the  simulations. 

[51]  03  loss  rates  are  generally  in  excess  of  100  ppbv/day 
within  the  rotational  lobes  of  the  vortex  at  this  time.  These 
regions  tend  to  be  more  isolated  due  to  relatively  weak  winds 
and  are  at  lower  latitudes  so  that  trajectories  within  the  lobes 
have  experienced  longer  periods  of  high  daylight  fractions 
than  trajectories  closer  to  the  pole  and  within  the  shear  zone 
(high  Q).  Peak  loss  rates  in  the  lobes  over  Greenland  and 
Northern  Scandinavia  are  as  large  as  140  ppbv/day.  The  air  in 
the  Scandinavian  lobe  was  sampled  by  the  ER2. 

[52]  The  comparison  between  the  observed  and  predicted 
CH4/O3  correlation  shows  that  the  LaRC  LCTM  does  a 
remarkable  job  in  characterizing  the  compact  nature  of  the 
CH4/O3  correlation  below  45  OK  (the  maximum  potential 
temperature  sampled  by  the  ER2  on  26/02/00  was  456K). 
The  measurements  (black)  overlay  the  model  points  with 
both  showing  O3  mixing  ratios  of  approximately  2.0  ppmv 
for  CH4  less  than  1.0  ppmv.  The  model  shows  lower  CH4 
mixing  ratios  at  450K  associated  with  the  region  of  strong 
descent  over  the  pole  that  was  not  sampled  by  the  ER2.  The 
lowest  observed  CH4  was  encountered  near  45 6K  during 
the  ER2’s  final  climb  before  descent  to  Kiruna,  Sweden. 
Consequently,  the  ER2  sampled  air  which  was  experienc¬ 
ing  relatively  high  loss  rates  and  was  well  removed  from 
the  region  of  strong  diabatic  descent  and  enhanced  diffu¬ 
sive  tendencies.  Above  450K  the  model  predicts  a  wide 
range  in  both  03  and  CH4  for  a  given  potential  temper¬ 
ature  bin. 

4.2.  Case  Study:  05/03/00 

[53]  The  second  case  study  considers  the  synoptic  sit¬ 
uation  during  the  ER2  flight  of  05/03/00  that  occurred  at  the 
beginning  of  the  large  predicted  03  loss  rates  (see  Figure  7). 
Figure  12  shows  synoptically  mapped  ensemble  predictions 
of  condensed-phase  HN03  (NAT),  Clx/Cly,  and  03  net 
photochemical  loss  rates  between  425K  and  475K.  A 
comparison  between  the  observed  and  predicted  correlation 
between  CH4  and  03  within  the  vortex  is  also  shown.  Only 


those  trajectories  with  Clx/Cly  ratios  larger  than  0.1  are 
shown  to  highlight  the  filament  of  processed  air  over  the 
Maritime  Provinces  of  Canada.  Streamline  contours  are 
shown  along  with  the  NAT  distribution  to  indicate  the 
counter-clockwise  circulation  within  the  vortex. 

[54]  The  model  predicts  a  PSC  encounter  on  the  eastern 
leg  of  the  ER2  flight  track  (shown  in  bold  black),  with 
condensed-phase  HN03  concentrations  reaching  6.0  ppbv 
within  the  core.  In  situ  NOy  and  particle  observations 
confirm  the  existence  of  this  PSC.  Photochemical  losses 
are  much  more  uniform  than  during  the  02/26/00  case  study 
since  much  of  the  vortex  is  outside  of  polar  night.  Photo¬ 
chemical  loss  rates  range  from  near  zero  over  Greenland 
(where  the  trajectories  have  just  come  over  the  pole  and 
have  spent  time  in  polar  night)  to  near  120  ppbv/day  or 
approximately  5%/day  along  the  eastern  leg  of  the  ER2 
flight  (in  the  region  of  the  synoptic  PSC).  Chlorine  activa¬ 
tion  is  very  high  (Clx/Cly  approaching  1.0)  within  the  PSC, 
while  just  to  the  southeast  of  this  region  Clx/Cly  ratios  fall  to 
near  0.6.  This  wide  range  of  activation  levels  is  reflected  in 
the  broad  modeled  Clx/Cly  frequency  distribution  during 
this  time  period  (Figure  9). 

[55]  The  variation  in  activation  levels  arises  as  trajectories 
orbit  around  the  vortex,  pass  through  the  synoptic  scale 
NAT  cloud  on  the  Siberian  side  of  the  vortex  where  they 
become  activated,  and  then  slowly  deactivate  as  they  move 
along  the  warmer  (near  2 1 5K)  Canadian  side  of  the  vortex. 
Trajectories  reach  their  lowest  level  of  activation  just  prior 
to  entering  the  western  edge  of  the  NAT  cloud,  which 
happens  to  be  south  of  the  ER2  flight  track  on  this  day. 
The  rapid  chlorine  activation  within  the  PSC  most  likely 
occurs  via  the  heterogeneous  reaction  C10N02  +  HC1  — > 
Cl2  +  HNO3,  since  there  are  ample  chlorine  reservoir 
species  available  just  upstream.  Highly  activated  trajectories 
are  also  found  near  the  center  of  the  vortex.  These  trajecto¬ 
ries  most  likely  passed  through  the  synoptic  PSC  on  some 
earlier  day,  or  were  activated  due  to  heterogeneous  reactions 
on  sulfate  aerosols  (predicted  vortex  surface  area  densities 
at  this  time  are  near  0.5  [im/cm 3)  and  have  remained 
activated  as  they  stagnate  in  a  region  of  relatively  cold 
temperatures,  limited  sunlight,  and  weak  circulation. 

[56]  There  is  a  filament  of  vortex  air  in  middle  latitudes 
extending  from  central  Siberia  to  the  west  coast  of  Green¬ 
land.  Examination  of  daily  maps  from  the  simulation  shows 
that  this  filament  is  a  remnant  of  an  earlier  planetary  wave 
breaking  event  that  occurred  on  29  February  1999,  and 
resulted  in  significant  export  of  processed  air  out  of  the 
vortex.  The  processed  air  rapidly  mixes  into  the  midlatitude 
background  due  to  strong  shear  deformation  on  the  equa¬ 
torial  flank  of  the  polar  night  jet.  However,  the  filament  still 
shows  03  photochemical  loss  rates  of  up  to  80  ppbv/day  and 
moderately  high  Clx/Cly  ratios  (0.4-0. 6),  5  days  after  it 
formed.  Planetary  wave  breaking  events  were  relatively 
infrequent  during  the  1999/2000  winter  (in  addition  to  the 
wave  breaking  on  29  February,  there  was  an  earlier  event  on 
10  February  at  450K).  In  general,  the  sharp  demarcation  in 
the  photochemical  characteristics  of  the  ensemble  trajecto¬ 
ries  (high  chlorine  activation,  high  O3  loss  rates)  is  coinci¬ 
dent  with  the  kinematic  edge  of  the  Arctic  vortex.  The 
coincidence  between  the  kinematic  edge  and  the  region  of 
photochemical  processing  arises  in  spite  of  the  fact  that  the 
photochemistry  is  largely  independent  of  the  edge  defini- 
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Figure  12.  Synoptically  mapped  ensemble  predictions  of  condensed-phase  HN03  (NAT)  (ppbv),  Clx/ 
Cly,  and  03  loss  rates  between  425K  and  475K  on  5  March  2000  (03/05/00).  The  lower  right  panel  shows 
a  comparison  between  the  observed  (black  triangles)  and  predicted  (colored  crosses)  correlation  between 
CH4  and  03  within  the  vortex.  The  black  line  in  the  synoptic  maps  indicates  the  kinematic  edge  of  the 
Arctic  vortex  at  45 OK.  The  bold  black  line  indicates  the  ER2  flight  track. 


tion,  highlighting  the  utility  of  the  kinematic  edge  in 
defining  the  limit  of  chemical  processing. 

[57]  The  predicted  correlation  between  CH4  and  03  once 
again  does  a  very  good  job  of  reproducing  the  observed 
compact  correlation  at  45 OK  where  the  measurements  over¬ 
lay  the  modeled  03/CH4  correlation  (03  near  1 .7  ppmv  for 
CH4  mixing  ratios  below  1.0  ppmv).  The  model  predicts  a 
less  compact  correlation  above  45 OK  with  a  wide  range  of 
CH4  and  03  within  relatively  narrow  potential  temperature 
ranges.  This  variability  is  similar  to  the  high  variability 
found  in  early  December  (see  Figure  2)  and  late  February 
(see  Figure  9)  and  suggests  a  relatively  inhomogeneous 


vortex  just  above  the  altitudes  where  the  maximum  03  loss 
rates  are  found. 

[58]  This  case  study  suggests  that  the  disagreement 
between  the  observed  and  modeled  chlorine  activation 
(Figures  10  and  13)  may  arise  due  to  the  detailed  timing 
of  chlorine  activation  and  de-activation  as  the  air  parcels 
orbit  within  the  vortex.  We  cannot  comment  on  the  history 
of  the  air  parcels  observed  by  the  ER2.  But  the  slow 
deactivation  predicted  after  trajectories  encounter  the  syn¬ 
optic  scale  PSC  on  5  March  arise  because  there  is  only 
about  0.01  ppbv  of  modeled  NOx  available  within  the 
vortex  during  early  March  (consistent  with  the  observed 
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Figure  13.  The  vertical  distribution  of  the  accumulative  sources/sinks  in  vortex  03,  NOy,  and  H20 
between  1  December  1999  (12/01)  and  28  March  2000  (03/28)  from  the  LaRC  LCTM  simulation  is 
shown.  The  upper  panel  shows  the  accumulated  effects  of  photochemical  production-loss  (P-L)  and 
vertical  advection  on  the  vortex  ensemble  mean  03  distribution  (ppbv).  The  middle  panel  shows  the 
accumulated  effects  of  Sedimentation  (Sed.),  evaporation  (Evap.)  and  vertical  advection  on  the  vortex 
ensemble  mean  NOy  (ppbv).  The  lower  panel  shows  the  accumulated  effects  of  Sedimentation  (Sed.), 
evaporation  (Evap.)  and  vertical  advection  on  the  vortex  ensemble  mean  H20  (ppmv). 


NOx  mixing  ratios  within  the  vortex),  and  many  of  the 
trajectories  in  the  warmer  Canadian  side  of  the  vortex  are 
still  in  darkness.  Interestingly,  by  12  March,  (near  the  end  of 
the  largest  chlorine  catalyzed  ozone  loss  and  just  prior  to  the 
predicted  rapid  deactivation)  the  processes  which  maintain 
the  model  Clx/Cly  ratios  have  changed.  At  this  time,  PSCs 
are  no  longer  present  but  most  vortex  trajectories  at  45 OK 
have  temperatures  near  200K,  indicating  that  heterogeneous 
chlorine  activation  on  sulfate  aerosols  is  likely  to  play  a 
significant  role  in  extending  the  high  chlorine  activation  out 
to  mid-March  in  the  model  simulation.  This  shift  in  the 
physical  processes  responsible  for  maintaining  the  high 


modeled  Clx/Cly  ratios  over  the  period  of  one  week  illus¬ 
trates  the  complex  interactions  between  the  temperature, 
sunlight,  and  heterogeneous  processing  history  of  vortex  air 
parcels  during  early  March.  It  is  not  surprising  that  we  see 
significant  differences  in  the  rates  of  the  predicted  and 
observed  chlorine  deactivation  given  the  complexity  of 
the  vortex  evolution  at  this  time. 

5.  Vortex  Budgets 

[59]  Figure  13  shows  budgets  for  vortex  03,  NOy,  and 
H20  in  the  LaRC  LCTM  simulation.  The  vortex  budgets 
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were  obtained  by  integrating  the  daily,  ensemble  mean 
tendencies  during  the  period  between  12/01  and  03/28.  As 
shown  earlier,  the  integrated  photochemical  03  loss  reaches 
a  peak  of  2.2  ppmv  at  45 OK.  The  photochemical  03  losses 
are  largely  (but  not  completely)  balanced  by  descent  of  high 
03  mixing  ratios  from  aloft,  which  is  near  2.0  ppmv  between 
425K  and  525K.  Sedimentation  results  in  a  peak  integrated 
NOy  loss  of  nearly  10  ppbv  at  500K  and  evaporation  results 
in  approximately  2  ppbv  increase  in  NOy  at  450K.  The 
integrated  change  in  NOy  due  to  vertical  advection  peaks  at 
500K  with  an  accumulated  increase  of  10  ppbv.  The  model 
water  vapor  budget  is  dominated  by  vertical  advection, 
which  brings  down  1.7  ppmv  of  H20  at  475K.  Sedimenta¬ 
tion  results  in  a  peak  integrated  H20  loss  of  0.6  ppmv  at 
500K  that  is  associated  with  the  late  December,  early 
January  ice  PSCs.  This  dehydration  is  very  close  to  that 
seen  by  in  situ  balloon  measurements  [Schiller  et  al .,  2002]. 
The  0.6  ppmv  integrated  H20  loss  at  420K  is  spurious,  and 
is  a  result  of  rapid  “rainout”  (sedimentation  at  temperatures 
above  ice  frost  point)  of  artificially  high  initial  water  vapor 
mixing  ratios  during  the  March  PSC  event.  The  high  initial 
water  vapor  mixing  ratios  occurs  because  of  uncertainties  in 
the  satellite  water  vapor  retrievals  when  NAT  aerosols  are 
present.  The  net  contribution  due  to  H20  evaporation  is 
very  small  in  the  model.  The  vortex  03,  NOy,  and  H20 
budgets  are  nearly  identical  in  the  simulation  without 
diffusion. 

6.  Impact  of  Mountain-Wave  Activity 

[60]  Mountain  waves  (MWs)  are  small-scale  gravity  and 
inertial  gravity  waves  generated  by  flow  over  topography 
that  can  propagate  into  the  stratosphere  under  suitable 
conditions  and  leading  to  large  temperature  perturbations. 
If  the  temperatures  within  the  MW  drop  below  NAT  or  ice 
condensation  points  than  orographic  PSCs  can  form.  Oro- 
graphically  induced  PSCs  were  ubiquitous  over  the  Norwe¬ 
gian  Mountains  during  SOLVE.  Carslaw  et  al  [1999]  has 
shown  that  MWs  can  be  a  significant  source  of  synoptic 
scale  PSCs,  extending  several  thousand  kilometers  down¬ 
wind  of  the  MW  temperature  anomalies. 

[61]  The  LaRC  LCTM  simulation  provides  a  unique 
opportunity  to  assess  the  impact  of  the  MW  activity  on 
ozone  loss  within  the  Artie  vortex.  The  impact  of  temper¬ 
ature  perturbations  associated  with  MW  activity  on  Arctic 
ozone  loss  is  investigated  by  introducing  daily  three- 
dimensional  temperature  anomalies,  predicted  by  the 
NRL  Mountain  Wave  Forecast  Model  (MWFM  2.0)  [Eck- 
ermann  and  Preusse,  1999],  into  the  isentropic  data  base 
used  for  the  LaRC  LCTM  simulations.  MWFM  2.0  is  a 
ray-based  extension  of  the  previous  version  of  MWFM  that 
considered  a  two-dimensional  hydrostatic  gravity  waves.  In 
the  new  version,  wave  group  velocities  are  calculated  using 
a  rotating  nonhydro  static  dispersion  relation.  Wave  ampli¬ 
tudes  are  calculated  using  conservation  of  vertical  flux  of 
wave  action,  with  rotational  and  nonhydro  static  terms 
retained.  Wave  breaking  is  accommodated  through  both 
convective  and  dynamical  instability  criteria.  MWFM  2.0 
uses  meteorological  model  forecasts  of  wind  and  temper¬ 
ature  to  generate  a  forecast  of  mountain-wave  activity.  The 
key  to  producing  this  forecast  is  a  detailed  representation 
of  global  topography,  which  includes  estimates  of  the 


locations,  widths,  heights  and  orientations  of  nearly  all 
mountain  ridges  on  Earth  with  scales  ranging  from  a 
hundred  kilometers  down  to  several  kilometers.  The 
MWFM  ridge  database  is  compiled  from  5x5'  digital 
elevation  data  using  objective  techniques  to  identify  sig¬ 
nificant  features.  This  approach  does  a  reasonable  job  of 
forecasting  mountain-wave  amplitudes  in  the  upper  tropo¬ 
sphere  and  lower  stratosphere  [Bacmeister  et  al ,  1 994]  and 
was  used  to  forecast  mountain-wave  effects  operationally 
for  flight  planning  during  SOLVE  (Eckermann  et  al., 
unpublished  manuscript,  2001). 

[62]  The  MWFM  provides  grid  point  estimates  of  the 
mean  amplitude  and  standard  deviation  of  the  small-scale 
temperature  fluctuations  using  the  Goddard  lxl  degree 
GEOS-2  Data  Assimilation  System  (DAS)  [Schubert  et 
al ,  1993],  winds  and  temperatures  at  standard  pressure 
levels  between  100  and  10  mb.  The  MW  statistics  are 
interpolated  to  the  LaRC  LCTM  isentropic  data  base  and 
then  daily  distributions  of  synthetic  MW  temperature  per¬ 
turbations  are  reconstructed.  At  each  grid  point  the  recon¬ 
struction  first  obtains  a  random  sign  of  the  MW  temperature 
perturbation,  which  multiplies  the  MWFM  mean  amplitude, 
and  then  adds  a  temperature  increment  obtained  from  a 
random  Gaussian  distribution  with  a  mean  of  zero  and  a 
standard  deviation  determined  by  the  MWFM  statistics. 
This  reconstruction  results  in  daily,  gridded  MW  temper¬ 
ature  perturbations  that  have  the  same  three-dimensional 
distribution  of  means  and  standard  deviations  as  the  MWFM 
predictions,  when  averaged  over  a  number  of  independent 
realizations  (days).  However,  each  large-scale  gridbox  has 
either  a  negative  or  positive  MW  temperature  perturbation 
on  a  given  day.  In  reality  each  trajectory  would  experience 
both  positive  and  negative  MW  temperature  perturbations 
as  it  transects  the  assimilated  grid  box.  This  means  that  the 
randomly  signed  MW  temperatures  will  underestimate  the 
true  MW  impact  since  50%  of  the  grid  boxes  with  MW 
activity  will  have  only  positive  temperature  perturbations. 
To  assess  the  bias  introduced  by  the  assumed  MW  temper¬ 
ature  perturbation  we  considered  two  different  MW  scenar¬ 
ios,  one  with  MW  temperatures  as  defined  above  and  one 
where  the  MW  temperature  perturbation  is  always  negative. 
The  randomly  signed  MW  temperature  perturbation  pro¬ 
vides  a  conservative  lower  bound  on  the  actual  impact  and 
the  negative  MW  temperature  perturbation  provides  an 
absolute  upper  bound  for  assessing  the  impact  of  MW 
activity  on  Arctic  03  photochemistry. 

[63]  Figure  14  shows  maps  of  the  monthly  mean  ampli¬ 
tudes  of  the  reconstructed  MWFM  temperatures  at  45 OK, 
along  with  an  indication  of  the  percentage  of  time  a  given 
location  was  within  the  Arctic  vortex  during  each  month. 
Peak  monthly  mean  MW  amplitudes  are  near  3. OK,  con¬ 
sequently  MW  perturbations  will  have  the  largest  impact 
where  the  large-scale  temperatures  are  already  relatively 
close  to  the  NAT  frost  point  (i.e.,  within  the  Arctic  vortex). 
Most  of  the  strongest  MW  activity  is  associated  with 
orography  over  North  America,  Tibet,  and  the  Alps,  which 
are  generally  located  outside  the  Arctic  vortex  and  within 
warmer  air  masses.  However,  regions  of  moderate  MW 
activity  are  predicted  within  the  Arctic  vortex  over  the 
southeastern  tip  of  Greenland,  Iceland,  Scandinavia,  Siberia, 
and  Alaska  with  monthly  mean  amplitudes  of  1.0-1. 5K.  Of 
these  regions,  Scandinavia  shows  the  largest  monthly  mean 
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Figure  14.  Maps  of  the  monthly  mean  amplitudes  of  the  reconstructed  MWFM  orographic  gravity- 
wave  (O-GW)  temperatures  (K)  at  450K  (shaded).  Contours  indicate  the  percentage  of  time  a  given 
location  was  within  the  Arctic  vortex  during  each  month. 


predicted  MW  temperature  perturbations  and  the  most 
continuous  exposure  to  the  large-scale  cold  pool  within 
the  Arctic  vortex.  In  fact,  during  each  month,  at  least  40% 
of  the  Scandinavian  MW  activity  was  within  the  Arctic 
vortex  at  45 OK  during  the  1999/2000  winter. 

[64]  The  MW  impact  is  determined  by  computing  the 
differences  between  the  accumulated  vortex  ensemble 
means  for  simulations  with  and  without  MW  temperature 
perturbations.  Figure  15  shows  the  impact  of  MW  temper¬ 
ature  perturbations  on  the  vortex  03,  NOy,  and  H20  budgets 
during  the  period  from  1  December  through  28  March  2000. 
The  left  panels  show  the  results  with  the  randomly  signed 
MW  perturbation  and  the  right  panels  show  the  results  with 
the  negative  MW  perturbation.  The  random  MW  perturba¬ 


tion  shows  small  (less  than  20  ppbv  or  2%)  increases  in 
photochemical  03  loss,  small  (less  than  0.15  ppbv  or  2%) 
increases  in  NOy  sedimentation,  and  moderate  (0.03  ppmv 
or  7%)  increases  in  H20  sedimentation  between  450  and 
550K.  Smaller  absolute  changes  in  NOy  and  H20  evapo¬ 
ration  occur  below  the  potential  temperature  of  the  peak 
changes  in  sedimentation.  The  negative  MW  perturbations 
lead  to  a  100  ppbv  (6%)  change  in  the  accumulated  03 
photochemical  loss  within  the  vortex  at  475K,  a  0.8  ppbv 
(8%)  change  in  the  accumulated  NOy  sedimentation  at 
500K,  and  a  0.15  ppmv  (18%)  change  in  the  accumulated 
H20  sedimentation  at  450K.  Consequently,  the  LaRC 
LCTM  simulations  indicate  that  the  impact  of  MW  on  the 
accumulated  photochemical  03  loss  within  the  Arctic  vortex 
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Figure  15.  Impact  of  mountain- wave  (MW)  temperature  perturbations  on  the  vortex  03,  NOy,  and  H20 
budgets  during  the  period  from  1  December  through  28  March  2000.  The  left  panels  show  the  results 
with  the  randomly  signed  MW  perturbation  and  the  right  panels  show  the  results  with  the  negative  MW 
perturbation. 


during  the  1999/2000  winter  is  a  2%  to  6%  increase  in  the 
accumulated  loss,  with  similar  impacts  on  the  accumulated 
NOy  sedimentation  and  7%  to  18%  increases  in  the  accu¬ 
mulated  H20  sedimentation. 

7.  Conclusions 

[65]  The  LaRC  LCTM  has  been  used  in  conjunction  with 
HALOE  and  POAM  III  satellite  observations  to  simulate 


the  large-scale  photochemical  evolution  of  the  Arctic  vortex 
during  the  1999/2000  winter.  Ensemble  averaging  of  the 
vortex  trajectories  allows  us  to  investigate  the  altitude-time 
distribution  of  the  vortex  evolution  while  synoptic  mapping 
of  the  trajectories  allows  us  to  conduct  focused  case  studies 
for  comparisons  with  in  situ  and  remote  data  obtained  by 
ER2,  DC8,  and  OMS  Balloon  payloads  during  SOLVE. 

[66]  The  LaRC  LCTM  ensemble  mean  vortex  predictions 
show  significant  denitrification  occurring  within  the  vortex 
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during  late  December  1999  and  early  January  2000  near 
500K  associated  with  extensive  PSCs  (primarily  con¬ 
densed-phase  HN03,  or  NAT)  within  the  vortex  at  these 
altitudes.  Large-scale  dehydration  due  to  sedimentation  of 
condensed  H20  occurred  briefly  during  this  same  period 
and  accounted  for  most  of  the  modeled  denitrification.  The 
model’s  ability  to  predict  the  combined  effects  of  diabatic 
descent,  diffusive  mixing  and  bulk  PSC  processes  appears 
to  be  reasonable  since  the  predictions  of  the  extent  of 
denitrification  within  the  vortex  are  relatively  good.  How¬ 
ever,  the  model  tends  to  underestimate  the  denitrification 
for  the  majority  of  the  vortex  by  approximately  10%,  and 
does  not  show  signatures  of  unprocessed  air  that  are  seen  in 
the  observations.  The  lack  of  unprocessed  air  in  the  model 
simulations  is  traced  to  overestimates  in  the  diffusive 
mixing  within  the  vortex. 

[67]  Predicted  chlorine  activation  reaches  over  80%  of 
the  available  Cly  during  late  January,  declines  slightly 
during  February,  and  then  increases  again  during  March 
as  renewed  PSC  activity  occurs.  The  observed  Clx/Cly 
ratio  is  near  0.4  during  the  March  ER2  flights,  signifi¬ 
cantly  lower  than  predicted  by  the  model.  The  synoptic 
case  study  on  05/03  shows  that  the  spring  time  model 
chlorine  activation  biases  arise  due  to  subtle  differences  in 
the  rates  of  chlorine  activation  on  synoptic  PSCs  and 
down-stream  de-activation  as  the  trajectories  orbit  around 
the  Arctic  vortex.  The  high  levels  of  activation  predicted 
on  03/05  are  maintained  because  of  low  NOx  and  nearly 
continuous  darkness,  while  mid-March  activation  is  main¬ 
tained  via  heterogeneous  processing  on  sulfate  aerosols. 
However,  the  exact  extent  of  the  overestimate  is  difficult 
to  assess  since  there  is  an  unaccounted  for  deficit  in  the 
observed  chlorine  species  of  approximately  33%  relative  to 
the  AC  ATS  total  inorganic  chlorine  measurements. 
Increases  in  the  observed  Clx  would  reduce  the  observed 
deficit  and  would  also  reduce  the  model  Clx/Cly  bias, 
while  increases  in  C1N03  or  HC1  would  also  reduce  the 
deficit  but  increase  the  model  bias. 

[68]  The  predicted  net  ensemble  mean  photochemical  03 
loss  reaches  over  60  ppbv/day  during  the  first  two  weeks  of 
March.  The  predicted  peak  loss  rate  is  in  very  good  agree¬ 
ment  with  that  inferred  from  the  MATCH  campaign  (—61.6 
±4.8  ppbv/day)  during  the  same  period,  resulting  in  small 
differences  in  the  accumulated  03  photochemical  loss  (—2.2 
ppmv  predicted  versus  —2.7  ±  0.24  ppmv  observed).  The 
majority  (55%)  of  the  modeled  photochemical  loss  is  driven 
by  C10x  catalytic  cycles  with  BrOx  cycles  contributing  the 
remaining  45%.  Vertical  advection  nearly  balances  the 
photochemical  03  tendencies. 

[69]  The  impact  of  mountain-wave  temperature  perturba¬ 
tions  is  investigated  by  blending  the  NRL  MWFM  moun¬ 
tain-wave  temperature  perturbations  into  the  LaRC  LCTM 
temperature  fields.  The  results  of  sensitivity  studies  show 
that  the  impact  of  mountain-wave  temperature  perturbations 
on  Arctic  03  depletion  and  NOy  denitrification  range  from 
2-8%  during  the  1999/2000  winter.  Larger  MW  impacts 
would  be  expected  during  winters  without  such  extensive 
cold  temperatures. 

[70]  This  study  highlights  the  combined  effects  of  differ¬ 
ential  diabatic  descent  and  horizontal  diffusion  on  the 
distribution  of  trace  gases  within  the  vortex.  Such  coupling 
is  unique  to  the  Arctic  because  transient  midlatitude  dis¬ 


turbances  tend  to  distort  the  Arctic  vortex  more  than  its 
Southern  Hemisphere  counterpart,  thus  leading  to  lobe 
structures  and  shear  zones  with  enhanced  small-scale  mix¬ 
ing  efficiencies  within  the  vortex.  Maximum  diabatic 
descent  rates  are  often  coincident  with  the  enhanced 
mixing  regions  since  the  midlatitude  disturbances  which 
cause  the  vortex  distortions  also  transport  warm  air  pole- 
ward,  leading  to  enhanced  radiative  cooling  near  the  polar 
night  jet  core.  The  LaRC  LCTM  simulations  show  that 
diffusive  mixing  is  necessary  to  accurately  predict  the 
distribution  of  chemical  species  within  the  Arctic  vortex 
during  the  1999/2000  winter.  However,  diffusive  mixing 
does  not  significantly  impact  the  ensemble  mean  03  loss 
rates  within  the  vortex. 

[71]  The  results  of  the  LaRC  LCTM  simulation  during  the 
1999/2000  winter  indicate  that  large-scale  models  with 
current  chemistry  can  do  a  reasonable  quantitative  job  in 
predicting  the  observed  03  loss  and  extent  of  denitrification 
in  the  Arctic  vortex  with  bulk  parameterizations  of  the  PSC 
microphysics.  However,  significant  disagreements  remain. 
First,  the  modeled  loss  rates  arise  as  a  result  of  sustained 
periods  of  high  chlorine  activation  during  March,  when  the 
modeled  chlorine  activation  is  significantly  higher  than 
observed.  Second,  the  extent  of  the  observed  denitrification 
is  underestimated  in  the  model  simulation. 

[72]  Acknowledgments.  Support  for  this  work  was  provided  by  the 
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